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SUMMARY 


This  report  presents  a  vector  theory  for  the  backscattering  of  electromagnetic 
waves  from  a  random,  rough  surface.  The  basic  technique  used  is  that  employed  by 
Dr.  Adrian  K.  Fung  in  an  earlier  work.  The  surfaces  used  in  the  report  are  those  that 
are  generated  by  a  stationary  gaussian  random  process  as  opposed  to  surfaces  gener¬ 
ated  by  a  random  array  of  objects.  The  former  type  of  surface  is  assumed  to  simulate 
many  of  the  vegetation-free  sections  of  the  earth,  it  is  important  to  understand  the 
basic  characteristics  of  scattering  from  such  surfaces  for  the  purpose  of  aiding  military 
geographic  analysis  by  radar. 
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GLOSSARY  OF  SYMBOLS 


o°  backscatter  coefficient;  ratio  of  differential  radar  cross  section  to  differential 
surface  area. 

y  backscatter  coefficient;  ratio  of  differential  radar  cross  section  to  differential 
projected  area. 

0  angle  of  incidence. 

d'  local  angle  of  incidence. 

X  incident  wavelength. 

k  wavenumber  in  the  propagating  medium  above  the  surface, 

k'  wavenumber  in  the  propagating  medium  below  the  surface. 

p  (x,  y)  function  describing  the  total  composite  surface, 
s  (x,  y)  function  describing  the  slightly  rough  surface. 

Z  (x,  y)  function  describing  the  large  undulations. 

n,  unit  vector  in  the  direction  of  the  incident  wave. 

7  range  vector  from  the  origin  of  the  (x,  y,  z)  coordinate  system  to  the  surface 
point. 

n  unit  vector  normal  to  the  total  surface,  p(x,  y). 

P  the  point  where  the  scattered  field  is  computed. 

R,  unit  vector  in  the  direction  of  P  from  the  origin  of  the  coordinate  system. 

Oj  unit  vector  in  the  direction  of  P  from  a  surface  point, 

to  angular  frequency  of  the  wave, 
p  permeability  of  the  medium  below  the  surface. 
e  permittivity  of  the  medium  below  the  surface. 

17  intrinsic  impedance  of  the  medium  above  the  surface. 

E  total  electric  field  on  the  surface. 


II  total  magnetic  field  on  the  surface. 

Et  scattered  field  at  point  P  and  in  the  direction  defined  by  the  unit  vector,  n^. 
R  the  distance  from  the  origin  to  the  point  P. 

-H 

E,  locally  incident  field. 

R|  Fresnel  reflection  coefficient  for  a  local  horizontally  polarized  plane  wave. 


GLOSSARY  OF  SYMBOLS  (Cont’d) 


C  (r) 

Ci(0 


Fresnel  reflection  coefficient  for  a  local  vertically  polarized  plane  wave. 
Fourier  transform  of  the  slightly  rough  surface  function  s  (x,  y). 
autocorrelation  coefficient  for  the  large  undulations, 
autocorrelation  coefficient  for  the  slightly  rough  surface. 


BACKSCATTERING  OF  ELECTROMAGNETIC  WAVES 
FROM  A  SURFACE  COMPOSED  OF  TWO  TYPES 
OF  SURFACE  ROUGHNESS 


I.  INTRODUCTION 

1.  Purpose,  The  purpose  of  this  report  is  to  present  a  general  vector  theory  for  the 
backscattering  of  electromagnetic  waves  from  random,  rough  surfaces.  The  surfaces 
under  consideration  will  be  those  generated  by  a  stationary,  random  process  as  opposed 
to  those  generated  by  a  random  array  of  objects. 

2.  Background.  At  the  present  time,  the  method  of  extracting  military  geographic 
information  about  terrain  features  from  radar  is  almost  solely  dependent  upon  the  qual¬ 
itative  analysis  of  radar  imagery  by  a  photointerpreter.  Radar  should  be  a  quantitative 
tool  in  view  of  the  controlled  frequency,  look  angle,  and  self-contained  illumination 
system.  In  order  to  improve  upon  the  present  process  and  to  make  it  more  quantita¬ 
tive,  an  understanding  of  the  process  of  electromagnetic  wave  interaction  with  and 
scattering  from  natural  surfaces  must  be  obtained.  A  general  solution  to  the  problem 
of  plane-wave  scattering  from  an  arbitrary  rough  surface  is  still  lacking.  However,  re¬ 
cent  advances  in  vector  approximate  solutions  have  yielded  considerable  insight  into 
the  scattering  phenomena.  Radar  scattering  theories  are  basically  separated  into  two 
broad  classes:  those  which  deal  with  surfaces  generated  by  stationary,  random  proc¬ 
esses;  and  those  which  deal  with  surfaces  generated  by  a  random  array  of  objects.  This 
report  will  deal  exclusively  with  stationary,  random  processes  because  they  are,  in  gen¬ 
eral,  more  realistic  of  natural  terrain  than  the  random  array  of  objects.  The  parameters 
that  affect  the  scattering  of  an  electromagnetic  wave  from  a  rough  surface  must  be  iso¬ 
lated  and  molded  into  a  vector  theory  which  will  allow  the  calculation  of  average  re¬ 
turn  power.  The  average  return  power  can  then  be  used  to  calculate  the  backscatter  co¬ 
efficients  o°  and  y  .  The  backscatter  coefficient  o°  can  be  defined  as  the  ratio  of 
differential  radar  cross  section  to  differential  surface  area,  while  y  can  be  defined  as 
the  ratio  of  differential  radar  cross  section  to  differential  projected  area.  A  simple  re¬ 
lation  exists  between  y  and  o°  which  is: 

o°  =  y  cos  0 

where  8  is  the  angle  of  incidence.  The  two  backscatter  coefficients  o°  and  y  are 
functions  only  of  terrain  properties  and  are  not  functions  of  radar  parameters.  One  of 
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the  most  difficult  problems  associated  with  rough-surface  scatter  is  the  question  of  de¬ 
polarization.  If  a  plane  wave  is  incident  onto  a  rough  surface  with,  say,  horizontal  po¬ 
larization,  it  is  possible  to  receive  not  only  a  horizontal  component  but  also  a  vertical 
component.  This  is  due  to  the  fact  that  the  rough  surface  will  scatter  some  of  the  ori¬ 
ginal  incident  field  into  a  different  polarization,  i.e.,  depolarization  takes  place.  This 
very  brief  initial  statement  on  depolarization  is  sufficient  to  permit  a  look  into  the 
fundamental  parameters  that  affect  the  scattering  of  an  electromagnetic  wave  from  a 
rough  surface.  Five  basic  parameters  listed  and  discussed  by  Cosgriff,  Peake,  and  Taylor1 
are: 


a.  Surface  Roughness 

b.  Angle  of  Incidence 

c.  Polarization 

d.  Complex  Dielectric  Constant 

e.  Frequency 

The  surface-roughness  parameters  would  include  such  things  as  the  probability  distri¬ 
bution  of  the  surface  heights  about  a  mean  plane  and  the  surface  autocorrelation  func¬ 
tion.  The  surface-roughness  parameters  are  the  predominant  factors  in  determining  the 
scattered  field.  A  qualitative  understanding  of  surface  roughness  can  be  obtained  through 
the  Rayleigh  criterion  which  relates  the  “roughness”  of  a  surface  to  wavelength  and  angle 
of  incidence.  To  this  date,  no  thorough  quantitative  measure  of  surface  roughness  has 
been  developed  which  may  explain  the  fact  that  a  general,  exact  solution  to  the  problem 
of  wave  scattering  from  rough  surfaces  is  still  lacking. 

The  angle  of  incidence  refers  to  the  angle  between  the  direction  of  the  incident  wave 
and  the  vertical.  The  plane  of  incidence  is  then  the  plane  containing  the  wave-propaga¬ 
tion  direction  and  the  vertical.  In  Fig.  1,  the  arrow  located  in  the  xz-planc  represents 
the  direction  of  the  incident  wave  with  angle  of  incidence,  0 ,  and  the  plane  of  incidence 
is  the  xz-plane.  At  the  far  ranges,  the  angle  of  incidence  will  be  large  which  usually  re¬ 
sults  in  difficult  problems  with  shadowing.  At  near  ranges,  the  angle  of  incidence  will 
be  small.  The  angle  of  incidence  at  a  point  on  a  surface  can  also  be  defined  as  the  angle 
between  the  direction  of  the  incident  wave  and  the  normal  to  the  surface.  In  this  case, 
the  angle  of  incidence  is  usually  referred  to  as  the  local  angle  of  incidence. 

The  polarization  of  an  electromagnetic  wave  is  defined  by  the  direction  that  the  elec¬ 
tric  field  vector  takes.  If  this  direction  is  constant  in  time,  then,  the  wave  is  said  to  be 
linearly  polarized.  This  report  will  consider  two  types  of  linear  polarization  in  particular; 


1  Coagriff,  Peake,  Taylor.  Terrain  Scattering  Properties  for  Sensor  System  Design  (T strain  Handbook  II),  Hie  Ohio 
Slate  Univenity,  Engineering  Experiment  Station  Bulletin  No.  181. 
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Fig.  1.  The  arrow  located  in  the  xz-planc  represents  the  direction  of  the  incident  wave 
with  angle  of  incidence,  0 ,  and  the  plane  of  incidence  is  the  xz-plane. 


horizontal  and  vertical.  Horizontal  polarization  occurs  when  the  electric  field  vector  is 
perpendicular  to  the  plane  of  incidence  and  the  magnetic  field  vector  lies  in  the  plane  of 
incidence.  Vertical  polarization  occurs  when  the  magnetic  field  vector  is  perpendicular 
to  the  plane  of  incidence  and  the  electric  field  vector  lies  in  the  plane  of  incidence.  Po¬ 
larization  effects  come  into  being  nut  only  because  of  the  direction  of  the  incident  field 
vector  but  because  of  the  depolarization  properties  of  the  surface.  The  complex  dielec¬ 
tric  constant  is  a  function  of  the  surface  electrical  properties  (permittivity  and  conduc¬ 
tivity)  and  the  frequency  of  the  incident  wave. 

Two  basic  techniques  have  evolved  for  working  the  problem  of  electromagnetic  wave 
scattering  from  a  surface  generated  by  a  stationary  random  process.  This  type  of  sur¬ 
face  is  one  in  which  the  elevations  or  heights  of  the  surface  above  a  mean  plane  are  dis¬ 
tributed  with  a  certain  probability  density  distribution  function.  The  distribution  func¬ 
tion  usually  assumed  is  gaussian  in  order  to  help  simplify  the  calculations.  The  tangent- 
plane  method  of  computing  wave  scattering  problems  consists  of  approximating  the  local 
fields  at  a  point  on  the  surface  by  the  fields  which  would  be  present  on  a  plane  tangent 
to  the  surface  at  the  desired  point.  The  calculated  surface  fields  are  then  placed  in  the 
Helmholtz  integral  for  the  scalar  case  or  the  Stratton-Chu  integral  for  the  vector  case. 
The  scattered  field  in  any  particular  direction  can  then  be  computed.  In  order  for  the 
tangent-plane  approximation  of  local  surface  fields  to  be  correct,  the  radii  of  curvature 
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of  the  surface  at  all  points  must  be  much  greater  than  the  incident  wavelength.  An  ex¬ 
act  geometrical  relationship  which  gives  the  criterion  for  the  tangent-plane  method  was 
developed  by  Brekhovskikn  and  can  be  stated  as  follows: 

4jrrc  cos0»  X 

where  r  is  the  smaller  of  two  radii  of  curvature  of  a  surface  at  a  point.  The  angle  O'  is 
the  local  angle  of  incidence,  and  X  is  the  incident  wavelength.  Therefore,  it  can  be 
easily  seen  that  the  tangent-plane  method  will  allow  only  the  use  of  smoothly  undulat¬ 
ing  surfaces  and  that  no  roughness  which  is  smaller  than  a  wavelength  can  be  tolerated. 

The  other  basic  method  of  solving  the  scattering  problem  is  known  as  the  small  per¬ 
turbation  method.  In  this  method,  the  Fourier  series  or  Fourier  transform  is  used  in  a 
perturbation  series  to  solve  for  the  unknown  fields  above  and  below  the  rough  surface. 
The  exact  boundary  conditions  are  used;  and,  theoretically,  this  method  would  yield  a 
good  solution  to  the  scattered  field  from  any  surface  if  enough  terms  of  the  perturbation 
series  could  1m  computed.  In  actual  practice,  it  becomes  exceedingly  difficult  to  take 
much  more  than  the  first  term  of  the  series,  particularly  if  finite  conduction  on  the  sur¬ 
face  is  allowed.  Since  the  perturbation  series  effectively  represents  the  results  of  the 
surface  perturbations  about  a  flat  plane,  then  only  slightly  rou^i  surfaces  can  be  tol¬ 
erated  with  this  method.  A  quantitative  definition  of  a  slightly  rough  surface  could  be 
as  follows: 

lks(x,y)|<  1 

where  k  is  the  wavenumber  in  the  propagating  medium  above  the  surface  and  s(x,y) 
is  a  surface  function  representing  elevations  from  a  mean  plane  to  the  surface. 

The  tangent  plane  method  has  been  shown  to  yield  good  results  for  the  average  back- 
scattered  power  with  the  like-polarized  component,  at  least  for  small  angles  of  incidence 
(0s -20°)  when  a  gaussian  autocorrelation  function  is  used.  When  an  exponential  corre¬ 
lation  function  is  used,  better  results  are  obtained  for  the  like-polarized  term.  It  can 
easily  be  shown,  however,  that  a  pure  exponential  correlation  function  leads  to  an  imagi¬ 
nary  surface  slope  distribution  so  that  results  using  this  correlation  function  are  very  mis¬ 
leading.  The  depolarized  components  from  the  tangent-plane  method  depend  upon  the 
sum  of  the  two  Fresnel  reflection  coefficients  which  does  not  compare  well  with  experi¬ 
mental  results  in  general. 

Many  natural  surfaces  can  be  assumed  to  be  made  of  large  undulations  with  small  per¬ 
turbations  superimposed.  Therefore,  a  correct  theory  on  the  backscattering  of  electro¬ 
magnetic  waves  from  rough  surfaces  should  consider  both  types  of  surface  roughness. 

2  Beckmann  wdSpiiHchina,  Tht  Scattering  of  Fire  tromagnt  tie  BWi  from  Rough  Surf  act  t,  Perga  mon  Pieae  1963, 
p.  29. 
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A  recent  paper  by  Fung  and  Haiaa-Lien  Chan1  has  done  thia  for  an  incident  wave  with 
horizontal  polarization  by  using  the  small  perturbation  method  up  to  the  first-order 
terms  as  an  approximation  of  the  fields  oo  the  surface.  When  the  zero-order  terms  of 
the  small  perturbation  solution  are  used  as  estimates  of  the  surface  fields  and  placed  in 
the  Stratton-Chu  integral,  the  tangent-plane  solution  results.  The  paper '  y  Fung  and 
Chan  did  not  consider  the  effects  of  along-track  surface  slopes  or  the  effects  of  higher 
order  slope  terms  in  the  power  calculation.  Our  purpose  here  was  to  apply  Fung’s  ideas 
such  that  the  effects  of  along-track  slopes  and  higher  order  slope  terms  are  considered  in  the 
horizontal  polarization  solution  and  to  perform  the  entire  solution  for  an  incident  wave 
with  vertical  polarization.  For  the  vertical  polarization  solution,  it  was  necessary,  first, 
to  solve  the  problem  of  a  vertically  polarized  wave  being  scattered  from  a  slightly  rough, 
dielectric  surface.  This  derivation  is  performed  in  Appendix  A.  In  the  case  of  horizontal 
polarization,  use  was  made  of  Rice’s  result3 4 5  as  derived  by  Fung.1 

In  practically  all  problems  concerning  the  scattering  of  electromagnetic  waves  from 
rough  surfaces,  certain  assumptions  must  be  made  in  order  to  make  the  solution  analyti¬ 
cally  manageable.  The  present  solution  is  no  exception.  Some  of  the  more  basic  assump¬ 
tions  which  were  made  are: 

( 1 )  Shadowing  effects  are  neglected. 

(2)  Only  the  far  field  is  calculated. 

(3)  Multiple  scattering  is  neglected. 

(4)  The  density  of  the  scattering  elements  is  not  considered. 

(5)  The  treatment  is  restricted  to  surfaces  that  arc  generated  by  stationary 
random  processes  and  which  yield  the  elevation  values  distributed  with  a  two- 
dimensional  gaussian  density  distribution  function. 

In  addition  to  the  above  assumptions,  certain  mathematical  assumptions  are  made 
which  allow  an  analytical  expression  to  be  obtained.  The  main  assumption  of  the 
tangent-plane  method  in  computing  the  fields  at  the  surface  was  not  used. 

The  derivations  will  be  preceded  by  a  discussion  of  the  geometry  of  the  problem  to  be 
solved.  Following  this,  the  Stratton-Chu  integral  will  be  introduced  and  modified  for 
future  calculations.  Next,  the  entire  scattering  problem  is  derived  for  the  case  of  a  hori¬ 
zontally  polarized  incident  plane  wave.  The  final  result  is  the  scattering  coefficient  o° 


3  A.  K.  Fung  and  Hsiae-Lien  Chan,  "  Baekseattering  of  Waves  by  Composite  Rough  Surfaces,”  IFF.F.  Transaction*  on 
Antenna  and  Propagation,  September  1969. 

4S.  0.  Rice,  “Reflection  of  Electromagnetic  Waves  from  SUghtiy  Rough  Surfaces,”  Communication!  on  Pare  and 
Applied  Mathematics,  VoL  4, 1951. 

5  A.  K.  Fung,  "Mechanisms  of  Polarised  and  Depolarised  Scattering  from  a  Rough  Dielectric  Surface,”  Joacnal  of 
the  Franklin  Institute,  VoL  2S5,  No.  2,  February  1968. 
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for  both  the  like  and  depolarized  terms.  The  same  problem  is  then  solved  for  the  case 
of  a  vertically  polarized  incident  plane  wave.  The  final  results  appear  as  graphs  of  o° 
versus  incident  angle  for  various  types  of  soils,  moisture  contents,  surface-roughness 
parameters,  and  frequencies.  A  discussion  of  the  final  results  and  their  possible  meaning 
concludes  the  report.  Throughout  the  report,  the  rationalized  MKS  system  of  units  is 
used.  The  notation  used  in  this  report  will  be  the  same  as  that  employed  by  Fung  and 
Hsias-Lien  Chan6  in  their  composite  surface  theory  so  that  the  two  results  can  be 
compared. 


II.  ANALYSIS 

3.  Geometry  of  Scattering  Problem.  A  plane  electromagnetic  wave  with  a  harmonic 
time  dependence  of  exp(jut)  is  incident  onto  a  rough  surface  described  by  a  function 
p(x, y)  .  This  function  p(x,  y)  represents  the  total  elevation  from  a  mean  plane  to  the 
surface  at  an  arbitrary  point  (x,  y).  The  total  surface  shall  be  composed  of  the  sum  of 
large  undulations  and  a  slightly  rough  surface.  The  large  undulations  shall  be  repre¬ 
sented  by  the  general  function  Z.(x,  y)  and  shall  be  required  to  conform  to  the  tangent- 
plane  method.  The  slightly  rough  surface  shall  be  designated  by  s(x,  y)  and  will  be  re¬ 
quired  to  conform  to  the  requirements  of  the  small  perturbation  method.  The  relation 
among  these  surface  functions  is  then  p(x, y)  =  Z(x,  y)  +  s(x,  y).  It  will  be  assumed  that 
Z(x,  y)  and  s(x,y)  are  generated  by  independent  stationary  random  processes  with 
zero  means.  The  average  values  of  Z(x,  y)  and  s(x,  y)  are  then  the  xy-plane.  The  ge¬ 
ometry  of  the  scattering  problem  is  given  in  Fig.  2.  A  right-handed  rectangular  Cartesian 
coordinate  system  is  set  up  near  the  surface  such  that  the  xy-plane  is  the  mean  surface 
height  or  elevation.  The  unit  vector  nj  is  a  vector  in  the  direction  of  the  incident  wave. 
The  angle  0  is  the  incidence  angle.  The  vector  T  is  a  range  vector  from  the  origin  of 
the  coordinate  system  to  a  surface  point.  The  unit  vector  n  is  a  normal  to  the  total  sur¬ 
face  of  p(x,^.  The  field  point  P  is  where  the  scattered  field  is  to  be  computed.  The 
unit  vector  R(  points  from  the  origin  of  the  coordinate  system  to  the  field  point  P. 

The  unit  vector  is  in  the  direction  of  the  field  point  P  from  a  surface  point.  When 
the  point  P  is  placed  in  the  far  field,  it  can  be  seen  that  rf3  =  Kj .  For  the  case  of  back- 
scattering,  ifj  =  -rTj  .  The  medium  above  the  surface  is  assumed  to  be  free  space.  The 
medium  below  the  surface  is  a  dielectric  with  zero  conductivity.  This  is  a  good  approxi¬ 
mation  for  many  soils  especially  at  high  radar  frequencies.  Another  right-handed  rec¬ 
tangular  Cartesian  coordinate  system  (  x,  y,  z)  is  set  up  with  the  origin  sitting  on  the 
large  undulations  Z(x,  y)  such  that  the  xy-plane  is  tangent  to  Z(x,y)  at  an  arbitrary 


*  A.  K.  Fun|  and  lUas-Lkn  Chan. 
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Fig.  2.  Geometry  of  the  scattering  problem. 


point.  The  z  -axis  is  then  normal  to  the  large  undulations.  The  propagation  constant  in 
the  medium  z<p(x,y)  will  be  designated  k'  and  is  equal  to  w  J tie  .  The  param¬ 
eters  inside  the  radical  sign  represent  the  permeability  and  permittivity,  respectively,  of 
the  material  in  the  medium  below  the  surface. 


The  scattered  field  in  the  direction  defined  by  the  unit  vector  is  then  ob¬ 
tained  from  the  Stratton-Chu  integral  as  stated  by  Silver:7 

E$  =  Krfj  X  j" (n  x  E  -  x  (ft  x  It))  exp(jkn2  *T)ds . 


(1) 


E  and  H  are  the  total  electric  and  magnetic  fields  on  the  surface  and  17  is  the  intrinsic 
impedance  of  the  medium  above  the  surface.  The  vector  n  is  normal  to  tiie  total  surface. 


-Jk  R 


R  is  the  distance  from  the  origin  to  the  point  P.  The  parameter  k  is  the  wavenumber 
in  free  space  and  is  equal  to  2nl\  where  X  is  the  wavelength.  All  quantities  are  known 
in  equation  (l)  except  the  two  surface-field  vectors.  The  surface  fields  can  be  written 
in  terms  of  the  local  coordinate  system  (x,  y,  z)  as  follows: 


7  Samuel  Silver,  "Microwave  Antenna  Theory  and  IM|n”  MIT  Rad.  I*b.  Semi  12,  McGraw-Hill,  1947,  p.  161. 
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E  =  l?Ej  +?Ef  +tET 

if=*Hr  +fHy  +«Hj 


where  are  unit  vectors  in  the  local  coordinate  system  (x,  y,¥).  The  surface 

field  components  Ej,  Err,  Ej-  are  in  terms  of  the  local  coordinates.  The  first  problem 
which  must  be  solved  is  to  obtain  3f,  y,  and?  in  terms  of  T,  X  and  IT  which  are  unit 
vectors  in  the  (x,  y,  z)  system.  The  surface  functions  will  no  longer  be  written  with  their 
(x,  y)  arguments: 


i  =  <-?z„  - fz,  +  r> < i  ♦  z;  +  z; >'*  «  {-Tz.-fZj+io 

dZ  37 

where  Zx  =  and  Zy  =  jy  .  These  terms  are  representative  of  surface  slopes  in  x 

and  y.  For  small  slopes,  the  term  Zjj  +  Z*  is  to  be  considered  small  with  respect  to 
one.  Throughout  this  report,  terms  involving  Z*  and  Z*  will  be  considered  negligible 
with  respect  to  terms  involving  Zx  and  Zy  .  The  direction  that  the  y-axis  takes  is  de¬ 
fined  by  the  unit  vectors  rfj  and  ?: 

f  =  (i  x  x  rfj  | 


?  =  y  x? 

rfj  =  T  sin0  -  IT  cos  0  . 

By  defining  the  local  coordinate  system  (x,  y,  z)  in  the  above  manner,  a  local  plane  of 
incidence  has  been  established  which  is  defined  by  the  vectors  rf}  and  ? .  It  can  be 
seen  that  this,  local  plane  of  incidence  does  not,  in  general,  coincide  with  the  original 
plane  of  incidence  defined  by  the  unit  vectors  it  and  £.  The  two  planes  will  only  co¬ 
incide  when  Zy  is  equal  to  zero.  The  local  angle  of  incidence  (0)  is  defined  as  the 
angle  between  ?  and  -it,  ,  and  can  be  computed  as  follows: 

cos  O  ' =  ”  n  |  •  i  *  cos  0  +  Zx  sin  0  . 

The  unit  vectors  y  and  x  can  be  computed: 

z  ••  n,  =  1 1  Zy  cos  0  +  j  (sin0  -  Zx  cos0)  +  kZ.  sin  0  | 

|z  x  n,  I  *  sin  0  -  Zx  cos0 
y  iZyDocos0  +  j  +  KZvl)osin0 
x  *  i  -  J  Zy  D0  cos0  +  k  7h 
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where  Do  **  (sin  0  -  ZI  cos  0 )"'  and  a  Zy  term  has  been  dropped  in  D0  .  The  term 
n  x  E  can  be  calculated  in  terms  of  the  surface  field  components  EL ,  E_  EL  and  the 
unit  vectors  i,j,andk: 

n  x  E  =  i  (Ey.  ZyIr0cosi?  -  Ey  -  E,  sy  J  +  J [Ej  +  Ey  Zy  D0  cos®  +  E,  sx  J 
+  Mp^-^Ey  ]  . 

The  subscripts  on  the  variables  s  and  p  refer  to  partial  derivatives  with  respect  to  the 
subscripted  variable.  The  term  n  x  H  can  be  written  directly: 

n  x  ft  =  i  [Hj  Zy  Docos®  -  H_  +  j  [II*  +  Hy  Zy  Doco«0  +  H,  ^  | 


+  Mp.H  —  P.Hyl- 

* 

In  the  case  of  backscattering,  the  term  n}  x  (n  x  E)  is  calculated 


ttering,  the  term  n2  x  (n  x  E)  is  calculated  as  -  n,  x  (n  x  E): 
-  r?(  x  (n  x  fi)  =  -  (i  cos  0  +  k  sin  ® )(E^-  +  Ey  Zy  l)0  cos  0  +  Ey  sx ) 

+  j  (E-  (ZyDoco8l  0  +  py  sin  9)  -  Ey  (cosfl  ■»  px  sin 9) 

-  Ely  sy  cos  9  |  . 


-  fcj  sy  cos  9  |  . 

The  final  term  to  Ik-  computed  is  n2  x  (n2  x  (  n  x  II)):  »2  x  (n2  x  ( n  x  11))  =  -  j  (II-  + 

H-  Zy  Dq cos 6  +  llf  sx  |  -  (icos0  +  ksin0)(ll-  (Z  I>H  cos*  0  +py  sin0)  -II-  (cos 0  +  px  sinfl) 
-  Sy  cos  9 )  .  The  above  terms  can  be  placed  in  tne  Stratton-Chu  integral  to  obtain  an 
expression  for  the  backscattered  field  in  terms  of  the  surfaee  field  components  in  local 
coordinates: 

Ef  =  K  J" J"  |  j  |  EL  Zy  l)Q  cosJ  9  -  Ely  cos®  -  E^  sy  cos  9  +  f? II- 

*  nlly  Zy  I)0  COS  9  +  ijIU  \  +  sin®  (py  V*-  -  p%  Ely  )) 

«•  *• 


f  i  eosfl  +  k  sin  9  W  EL 


-  nil-  (Zy  D^cos*  9  +  py  sin®)  +  t?  ll_(cos  9  +  p%  sin  0) 

+  ijllj  8y  cos®  |  |  exp  (  -  jkn  t  •  f)  dy  dx  .  (2) 

The  limits  of  integration  are  over  the  illuminated  area.  The  only  unknowns  in  equation 
(2)  are  the  six  surfaee  field  terms  expressed  in  local  coordinates.  In  order  to  proceed 
further,  a  particular  incident  polarization  must  he  chosen  and  then  the  local  fields 
calculated. 


* 


4.  Horizontal  Polarization. 

a.  Like-Polarized  Term.  A  horizontally  polarized  plane  wave  with  unit  ampli¬ 
tude  is  incident  onto  the  statistically  rough  dielectric  (zero  conductivity)  surface  p(x,  y). 
It  shall  be  required  to  obtain  analytical  expressions  for  c/H  H ,  the  scattering  coefficient 
for  the  like-polarized  return,  and  for  0°  v  ,  the  scattering  coefficient  for  the  depolarized 
return.  In  order  to  do  this,  an  expression  must  first  be  obtained  for  EHH ,  the  like- 
polarized  backscattered  field.  From  this,  the  mean  power  density  <  EHH  E*H  >  must 
be  calculated.  Similarly,  equations  must  be  derived  for  EHV  and  <  EHV  E*v>  in  order 
to  calculate  0°  v  .  The  brackets  are  here  used  to  denote  the  fact  that  an  average  of  the 
quantity  inside  is  being  taken.  The  asterisk  is  used  to  indicate  the  complex  conjugate. 

0°H  can  then  be  calculated  from  the  equation 

„o  _  lim  4jt  R*  <  Ehh  > 

°«H  "  R«  “  - - 

Where  E,  is  the  value  of  the  incident  fieljl  without  regard  to  polarization.  The  same 
method  can  be  used  to  determine  0“  v .  The  incident  wave  for  horizontal  polarization 
can  be  written  as 

Ej  =  j  exp  |- jk(xsin0  -  zcos  0)  f  . 

The  time  harmonic  of  exp  (jut)  has  been  suppressed  and  will  not  be  carried  along  any 
further. 


in  order  to  determine  the  surface  fields  in  local  coordinates,  the  incident 
wave  must  first  be  written  in  terms  of  the  local  coordinates.  The  locally  incident  field 
E[  is  then 

=  |x  (x  •  j )  +  y  (y  •  j  )  +  z  (z  •  j  )j  exp  |-jk(x  sin0  -  Z  cos0)  f  . 

exp{-jk|(n,  •  x)T  +  (n,  •  y)y  +  (n,  •  z)z  |  • 

The  quantities  x  and  Z  represent  the  coordinates  of  the  origin  of  the  local  coordinate 
system.  The  original  polarization  of  the  wave  has  t>een  broken  up  into  thiee  compo¬ 
nents  in  the  local  coordinate  system.  The  above  dot  products  can  be  easily  calculated 
from  results  already  derived: 

x  -  j  =  -  Zy  Docos  0 

y  -j  =  1 

1  *  J  ly  ’ 
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Since  y  •  J  =  1,  and  this  ia  the  magnitude  of  the  original  tranamitted  field  vector,  the 
effectaofthe  x  •  j  and  the  z  •  j  components  are  neglected.  This  should  be  a  very 
good  approximation  for  small  Zy  which  "must  be  the  case  in  order  for  the  tangent-plane 
method  to  be  applicable  to  the  smoothly  undulating  surface.  The  quantities  rf,  •  x, 
n,  •  y  ,  and  nj  •  i  ate  easily  computed  from  the  previous  work: 

n,  •  x  =  sin  0  -  Zx  cos  0 

».  y  =  IZr  cos0  sin  0  -  ZyCosfl  sin0]/(sin0  -  Zxcos0)  =  O 
it,  -z  =  -  cos0  -  Zx  sin0 

The  local  angle  of  incidence  ( 0 )  can  be  related  to  angle  0  by  the  following  expressions: 
cos  0'  *  cos0  +  Zx  sin  0 
sinfl'  sin#  -  Zx  CO3  0 

Using  the  expressions  abo .  e  for  the  local  angle  of  incidence,  the  local  incident  field  can 
be  written  as 

E1,  *  y  exp  j- jk(xsin0  -  Zcos0)|  exp  |-jk(xsin0'- Ecos0')|  . 

The  problem  locally  is  njw  one  of  a  horizontally  polarized  plane  wave  incident  onto  a 
slightly  rough  dielectric  surface.  This  problem  has  been  solved  previously  by  Rice*  and 
Fung.9  The  fields  up  to  first  order  in  perturbation,  as  expressed  by  Fung  in  the  Fourier 
transform  notation,  will  be  the  form  used  here: 


V  iff 


ky  kx  Q'OTdkxdkS  EXF 


^exp(-jkx  sin0')£exp(jkz  cos 0)  +  R^exp(-jkz  cos0  )J 


DO  OO 

-  JJ  <k|  +  kx  k'? )  QTKP  dkx  dky  \  EXF 

.M  .AS  e 


*S.  0.  Kice. 

*  A.  K.  Fung,  “Mectuu.—-. 
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(3) 


Q-  _  jT(k,a  -  k* )  S  (k,  +  kainfl;^) 

"  2w  (kx  +  k'  )(k*  +k’  +  l^k;) 

T  =  1  +  Rl  . 

=  Fresnel  reflection  coefficient  for  a  local  horizontally  polarized  plane 
wave. 

EXP  =  exp  [-jk(xsin0  -  Z  cos  0)) 

EXP  =  expljk^x  +jkyy  -jk,* ) 

kx  ~  J ka  -  k|  -  kj  when  kJ  >  ka  +  kj 

kf  =  -  j^kjj  +  ka  -  k1  when  k|  +  ky  >  kJ 

k't  =  ^/k'1  -  k*  -  k*  when  k'a  >  k*  +  k* 

k't  =  -j  J  k]  +  ky  -  k'1  when  kjj  +  k*  >  k'J 

S(kjt ,  ky )  is  the  Fourier  transform  of  the  sli^itly  rough  surface  function 
s(x,  y)  with  k%  ,  ky  as  the  Fourier  variables: 


It  should  be  noticed  that  the  kx  argument  in  Sfk^ky)  has  been  replaced  by  kt  + 
k  sin  0’in  Q’.  The  angle  0  represents  the  local  angle  of  incidence  and  is  the  angle  be¬ 
tween  the  unit  vectors  z  and-n(  .  The  Fresnel  reflection  coefficient  is  in  terms  of 
the  local  coordinate  system  and  can  Ik:  written  as  follows: 

kcos0- ^/k^  kJ  sin*  0' 

1,1  =  kco,*'*  A,J  "-V  rin'T 
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In  place  of  cos®’  and  sin®',  the  following  expressions  are  used: 

cos®'  *»  cos®  +  Z,  sin® 
sin®’ m  sin®  -  L%  cos®  . 

The  present  form  of  Rj^  would  make  future  calculations  very  difficult  due  to  the  Lx 
in  the  radical  and  in  the  denominator.  It  would  be  advantageous  to  approximate 
with  the  first  two  terms  of  a  Taylor  series  expansion  about  l%  =  0  since  small  slopes 
are  assumed.  When  this  is  done,  R^  can  be  written  as 

Ri  ^  K  +  gz, 

where  R  = 

O 

g  =  -2k  Ro  sin®/  k'cos$ 
cos  $  =  -  k1  sin1 0  /k* 

Hie  angle  0  would  be  the  angle  of  refraction  if  the  surface  were  flat,  e.g.,  p(x,  y)  =  0. 
The  angle  4>  is  then  related  to  the  angle  of  incidence  0  by  Snell’s  Law.  The  magnetic 
fields  associated  with  the  electric  fields  given  by  equation  (3)  can  be  computed  from 
Maxwell’s  equation  written  in  the  local  coordinate  system: 

rjH-  jexp(-jkX8in0)|exp(jkz  cos®)  -  R^expf-jkz  cos®))  •  cos®' 


<ka  -  k'a )  Q'  EXP  dkE  dky 


EXP 


sin®'exp(-jkx  sin®)|exp(jkzco80')  +  RA  exp(-jkz  cos®‘)j 


Q  EXPdk.dk,}  EXP  . 
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At  the  surface,  where  the  fields  must  be  evaluated,  the  x  and  y  coordinates  are  zero 
and  J  **  s(x,y) .  The  term  Q'  can  also  be  written  as 


0  =  (1  +  R0  +  )  QS(kx  +  k  sin 9’,  ky ) 

with  Q  = 


j<k'a-k’) 

1  A  1.3 


2a(kl+k'f)(l4+lS|,  +  klk'i) 


Hie  fields  evaluated  on  the  surface  then  become 


= 


Ev- 


oo  or 

{(1  +  Ro>  //^,QSexp(-jMdMk  + 

-OO  .oo 

BZ,^/*kykxQsexp  (-jk^dl^  dk^EXP 

•oo  .oo  *  | 

=  |  cxp(jkscos0)  +  Rc  exp  (- jks  cos0)  +  gZx  exp(-jks  cos0) 

OO  OO 

-  0  +  R„)  ff  (ka  +kik'i)QSexp(.jkts)dki 


A 


-BZ.  +  Kkj  QS  exp(-jkts)dki  dky  |  EXP 

-OO  .OO 

E*  =  -  ^0  +  R*) yyikyk',QSexp(-jkls)dki(dky  + 

&*jf  kyk,QSexp(-jkis)  dk,  dkyJ  EXP 
nllj  =  |  cosfl  exp(jks  cos0)  +  Zxsin0  expfjks  cosff)  - Rc  cos0  exp(*jks  cos0) 
-  gZyco80  exp(-jks cos0)  -  Rq  Zy  sin0  exp(-jks  cos0) 


QSexpf-jkjSjdk^dky 


♦  gZ,  ff  [*s.<k.  -1- -  kr  k.  ]  QSexp<.jkI.>dJt„dk,}  EXP 

-  oo  _oo 

M  M  10 

nH,  =  I  (1+R„)  J  J  ^  *  (V,  -  k’i)OSrxp(-jki*)dkidki 

+  *//  ^<k.  -  U  OS«p<-jkl,)dki  Jky  J  EXP 

•ftl  -OO 

ijH.  =  |  sin0exp(jks cosfl)  -  Zx  cos#  exp(jks  cosfl)  +  Rosin0  exp(-jks cos®) 

-  Ro  Lx  cos#  exp(-jks  cosfl)  +  gZx  sin0  exp(-jks  cosfl) 

ff[K(K+K  +  k,k') 

+  (1+R°)  JJ  J - £ - —  QSexp(-jks)dkxdky 

.OP  -OO 

OO  OO 

C  C  k  fk*  +  k*  +  k  k  i  \ 

+  g\J J  - i-^QSexp(-jkt8)dk1|dky  }  EXP  .  (4) 

•oo  .00 

The  parameter  S(kx  +  k  sin  flC  ky )  has  been  written  without  its  arguments.  The  J  com¬ 
ponent  of  equation  (2)  for  Ef  together  with  the  above  equations  for  the  surface  fields 
can  be  used  to  calculate  EHH  ,  the  like-polarized  backscattered  field  for  a  horizontally 
polarized  incident  wave.  E1HH  in  terms  of  the  surface  field  components  is,  from  equa¬ 
tion  (2): 

Khh  =  K  L)o  COS1  6  -  CO80  -  Ej  8y  cosfl  +  Tilly 

+  Tilly  Zy  D0  cosfl  +  tjHt  hx  +  sin 6  py  Eg 

-  .infl  px  Ey  ]  exp(-jkn,  *r)dy  dx  . 
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When  the  surface  field  equations  (4)  are  placed  in  the  above  expression  for  EHH ,  the 
following  equation  results: 


Ehh=-2*oK 


J (costf  + 


goZ]|sin9j  exp(-jkn|  •  f)  EXPexp(-  jks  cos0)dydx 


_  j  2k  cos0 
-  k'  cos  <(> 


*hh  "  K  ffSf  {A  +  BZ*  +C8*  +  DZy  +Esy  }  SexP<-jk,  s)  EXP  exp(>jknt  -  r) 


dydxdk  dk 


f  k  (k2  +  k  k  )  k2  k  1 

where  A  =  (I  +  RJQ  cos0  (k2  +  k.k.)  +  (I  +  R„)  Q  |  ‘  \  ‘  -  +  -j^J 

,  fk,(k2 +k  k’>  k2k’  1 

B  =  g  cos0  (1^  +kik^)Q  +  g  [ - £ - +  J  Q 


+  (1  +  R.o)sin0(tfa  +kfkJQ 

f  k  (k2  +  k2  +  k  k’ )] 

C  =  (1+RJ  [ -*  « . - !_LjQ  +  (|+Ro)8in0  (k*  +  k,k,)Q 

I)  =  (1  +  Rq)  ky  kx  Q  cscffcoifti  +  csc0  cos0(I  +  Ro)  [  (k,  *  kj)  j  Q 

+  (1  +  Ro)sin0kyk)t  } 

E  =  ( I  +  R0 )  cos0  kykrQ  +  (l  +  Ro )  sin# kykxQ  . 

In  the  .ibove  calculations,  the  term  D  was  approximated  by  the  first  term  of  its  Taylor 
series  expansion  around  Zx  -  0.  The  reason  more  terms  of  the  series  are  not  taken  is 
that  higher  order  slope  terms  would  result  in  IHH  and  these  are  assumed  negligible. 

This  results  in  the  following  approximation  for  D  : 


D  m  esc  0  . 

O 
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This  approximation  will  be  poor  at  very  small  angles  of  incidence  but  should  be  accept¬ 
able  over  a  large  range  of  incidence  angles,  particularly  when  the  slopes  are  not  large. 
Hie  range  vector  r  is 

r  =  ix  +  Jy  +  kp(x,y)  =  ix  +  jy  +  k(Z  +  s). 


Using  the  above  expression  for  the  range  vector,  the  backscattered  field  EHH  becomes 


(a,  +  b,  7j%  )  exp  |-2jkx  sinfl  +  2jkZ  cos0  ]  dydx 


+ 


1 


Hll 


a,  =  2Ro  cos0  b,  =  2Rogo  sin0 
1HH  =  K  JJJJ |  A  +  BZ^  +  Cs^  +  DZy  +  Esy|  S  exp(jas)  exp(-2jkx  sii 


sin  0 


+  2jkZ  co80)dydxdkji  dk 


a  =  k  cos 0  -  ky  . 


It  is  interesting  to  see  that  EHH  consists  of  the  sum  of  two  basic  terms  the  first  of  which 
depends  solely  on  the  large  surface  undulations  Z(x,y)  to  which  the  tangent  plane 
method  applies  and  represents  the  quasi-specular  component  of  the  scattered  field.  The 
term  IHH  arises  because  of  the  small  surface  perturbations  s(x,y);  although  it  also  con¬ 
tains  terms  dependent  upon  the  large  surface  undulations.  In  calculating  EHH  E^  ,  it 
will  be  assumed  that  the  two  terms  are  independent  io  that  cross-product  terms  would 
produce  a  negligible  effect  upon  the  final  result: 


Ehh E™  =  K K"  *  *, b, \  *  *, b, K  +  b’, z. z. }  • 

exp  |-2jk  sin0(x  -  x  )  +  2jk  cos0  (Z  -  Z)|dydxdy'  dx' 


+  I 


HH 


1  • 
*HH 


In  order  to  compute  <EHH  E*H>  ,  a  two-dimensional  distribution  must  be  assumed  for 
both  Z(x,y)  and  s(x,y).  The  distributions  to  be  used  for  both  surfaces  will  be  gaussian 
with  zero  means.  The  variance  of  the  large  undulations  Z(x,y)  will  be  designated  o 1  . 
The  variance  of  the  slightly  rough  surface  s(x,y)  is  given  the  symbol  o*  .  With  a  two- 
dimensional  gaussian  distribution,  not  only  must  the  mean  and  variance  be  known  but 
the  normalized  autocorrelation  function  of  the  surfaces,  called  the  autocorrelation  coef¬ 
ficient,  must  be  given  a  particular  form.  The  symbol  for  the  autocorrelation  coefficient 
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for  the  large  undulations  will  be  C(r)  and  that  for  the  slightly  rough  surface,  C,  (r).  The 
form  of  C(r)  will  be  that  suggested  by  Fung: 10 

^  ’  CXP  {  L*(G  +  r/L)} 
r  =  A  *  x')1  +  (y  -  y)1 


x  -  x’ 


■  u  =  r  sin  a 


y  -  y  -  v  =  r  cos  a  . 

It  can  be  seen  that  C(r)  is  a  two-parameter  correlation  function,  and  when  G  is  zero  an 
exponential  function  results.  The  quantities  L  and  G  are  physically  significant  in  that 
they  combine  to  give  a  measure  of  the  distance  between  hills  or  valleys  on  the  surface. 
With  a  large  L  or  G,  the  distance  between  hills  will  be  large.  When  G  is  small,  C(r) 
will  be  essentially  exponential  except  near  the  origin,  and  this  will  result  in  a  physically 
realizable  surface  with  the  variance  of  the  surface  slopes  o*  being 

o  =  -o  C  (0)  =  ~J~  ■ 

L  u 

The  primes  represent  derivatives  with  respect  to  r  .  The  form  of  C,  (r)  will  be  assumed 
gaussian  and  can  be  written  as 

C,(r)  =  e"J',J  . 

The  parameter  I  is  the  correlation  distance.  The  following  averages  are  computed  in  Ap¬ 
pendix  B  and  are  needed  to  calculate  <EHH  EJH>  . 

<exp  (2jk  cos  0  (Z  -  Z))>  =  exp  |- Kf  (1  -  C)j 

<  Z|exp(2jk  cos0(Z  -  Z))>  =  <Zx  exp|2jk  cos0  (Z  -  Z)|> 

=  -  j2koJ  cos0  exp  |-  K,  ( 1  -C)| 

<Z* Z^exp  [2jk cos0  (Z  -  Z)|>  =  -  oJ  j  +  K,  (^J  ]  exp|-K((l -C)| 


10  A.  K.  Fung,  Correspondence  to  author,  November  9, 1970. 
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Kj  =  4k2  o2  cos2  6  . 


The  correlation  coefficient  has  been  written  without  its(u,v)  arguments.  Placing  the  aver* 
ages  in  the  integral  results  in  the  following: 


<EhhEh#h>  =  ™ 


lJ  ,[  *C 
-  b ,  <r  rr  + 
>  du 


'■m\) 


-4ja,b,  k <?  cos#  4— 
•  \  \  9u 


exp  (-2jku  sinS)  exp  [-Kj  (1  -C)J  dudv 


+  <^hh*h*h:>  • 


It  will  be  easier  to  evaluate  the  above  integral  if  the  coordinates  are  changed  from  the 
rectangular  (u,v)  to  the  polar  (r,a).  Realizing  that  the  Jacobian  is  equal  to  -  r,  the  ex¬ 
pression  for  <EHHE*H>  now  becomes: 


<Ehh  Ehh 


(§))}’ 


L  2  ,  8C 

ko*  coeJ  - —  sina 
o  T 


r  exp  (-2jkr  sin  a  sin0)  exp  |-K,  (1  -  C)]  drda 


+  <lK»li»>  • 


When  K(»  1,  a  very  rough  surface  results,  and  it  can  be  seen  that  the  main  contribu¬ 
tion  to  the  integral  in  r  comes  only  from  the  region  about  r  =  0  .  The  limits  on  r  can 
then  be  made  from  zero  to  infinity.  The  limits  on  the  integral  in  a  will  be  from  zero  to 
2>r  .  The  correlation  function  C(r)  can  be  approximated  by  taking  the  first  two  non¬ 
zero  terms  of  its  Taylor  scries  expansion  about  r  =  0 .  C(r)  then  becom<»« 


C(r)  * 


Taking  the  indicated  partial  derivatives  of  C(u,v)  with  respect  to  u  yields  the  follow¬ 
ing  expression  for  <  EHH  E  *H  >  : 
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<EHH  EHH 


If  M 

o  o  V 


-  16b*  o4  k*  cos*  0 1*  sin?  al 


lr rnnft  ^  °* 

kco‘'  "W  +  TcT 


•k  y 


^exp  (-2jkr  sin#  sina)  exp|'  j  r  dr  da 


In  order  to  perforin  the  integration  in  a,  the  following  integrals  are  needed  (see  Appen¬ 
dix  C): 


r 

I  exp  (-  2jkr  sind  sina)  da  =  2n  Jc(2kr  sin0) 


sina  exp(-2jkr  sintf  sina)  da  =  -j2ir  J,  (2kr  sin0) 


o 

/(  J  (2kr  sinff)  1 

sin  a  exp(-2jkr  sin0  sina)  da  =  2ir  ( — 2k"r"rin0 —  -  J2(2kr  sin0)  1  . 


The  J  functions  represent  Bessel  functions  of  the  first  kind: 


<EhhEhV  =  KK* 


:/{2”’ 


4irb2  o2 


16*  a  b  k  cob0  o*  rj  (2kr  sin  0) 
{  Jo(2krsin0)  +  - 1 - 

_ _  32sbj  a4k*  cos20r2  f  J|  (2kr  sin0) 

+  “Pn  -Jo(2krsin0)-  uj7<n  [  2krsin0 


-  J2(2krsin0)j  jexp[--|^]rdr  +  <1H„I*„> 

The  integration  in  r  can  be  carried  out  with  the  aid  of  the  following  integral: 


The  following  expression  is  then  obtained  for  <EHH  E*^  : 
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/  **l L* 

iE.„>  =  Kk.|-l_ 
L,Gb,1o,k,gin,fl  1 

~ — rxp 


G  8#oa  LlG«j  bj  k1  sin0co80 
"+ 


exp  (-kVG.inVKJ +<!„„!•„>  . 


The  quasi-specular  term  has  been  determined.  The  component  <IHH  IJH>  must  now  be 
calculated  to  complete  the  result.  From  previous  work,  an  expression  for  IHH1^H  can 
be  written 


‘hhIJh  =  KK 


•////////(  AA*  +  AB*^  i  AC*sy  +  AD*Z \ 


+  AE*s;  +  BA^  +  BB*Z>  t%  +  BC*Z^  +  BD*ZxZy  +  BE#Zx8y 
+  CA#^  +  CB*^  Z^  +  CC\ \  +  CD*8s Z'y  +  CE-s^  *' 

+  I)A*Z  +  DB*Z  t  +  DC*Z  s'  +  DD*Z  Z  +  DE*Z  s* 

i  *  “  y  •  y  y  yy 

+  EA*sy  +  EB*sy  tx  +  EC*sy  s’  +  ED*  ^  Zy  +  EE*^  s*y  |  . 

SS*  exp  (jas  -  ja*s’)  exp  2jk  sin0  (x  -  x)  +  2jk  cos0(Z  -  Z)|  . 
dy  dx  dydx’dk^  dky  dk,  dk’y  . 

In  order  to  calculate  <IHH  ,  it  is  necessary  to  use  the  following  averages  which  are 
computed  in  Appendix  B: 

<ZxZ’y  exp  |2jk  cos0(Z  -  Z)J>  =  <ZsZy  exp  1 2 jk  cosO  (Z  -  Z))> 


<SS*  exp(jas  -  ja*s’)>  *  2s6(k'i  -kx)6  (ky  -ky)ot5  Wfk^  +ksin0,ky) 


,|4o;  | a*  +  a#l  *  2  aa*!;,!^ 


<exp  [2jk  cos 0  (Z  -  Z)j>  =  exp  |-K,  (l  -  C)| 
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<Zt  exp  (2jk  coafl  (Z  -  %)]>  =  <%x  exp  (2jk  cosfl  (Z  -  Z)l> 

=  -  j2koJ  coaO  ^exp(-K,(l -C)) 
<Zy  exp[2jk  cos0  (Z  -  Z)]>  =  <Z'y  exp  [2jk  cosff  (Z  -  Z))> 

=  -  j 2 k  o1  cosfl^J  exp  [ -  K,  ( 1  -  C)) 

«  30, 

<^SS*exp(jas-ja*8')>%  -ja*o,  —  2n 6  (k',  -  kM )  6  (ky  -  ky ) 


W  (kx  +  k  sinfl ,  ky )  exp  ^  o2  (a2  +  a*1  •  2aa*C, 


dC  ' 

<8^8*  exp(ja»-  ja*s’)>  *  -  jao*  ~  2rr  fi  (k'^  -  km  )6  (ky-ky) 

•  W(kx  +  k  sinfl,  ky)exp^-^  o2  (a2  +  a*2  -  288*0,  )| 

4  ac, 

<8ySS*exp(jas  -  ja*s')>  -ja*o,  2»6(kx  -  ka)6  (k^,  -  ky) 

•  W(k<  +  k  sinfl,ky)exp|-  ~  o]  |«*  +  a*’  -  2  aa*C,  | 

<8*s,ss#exp(ja8'ja#8>  %  -2#a*  6(ky  -k^Stk'  -  k  )  W  (k  +  k  sin  0,  k  ) 


a*2  -  2  aa*C. 


<Zy  Zy  exp  |2jk  co80  (Z  -  Z))>  =  -  o 


•  {it1  +  ’>  (^r) }  ,xp{  t 

+  a*2  -  2aa*0,  |J 

,)l>  =  v{l7  +  *■(£)’}  "ipl 


2  i  2 

o  |  a 


exp  |-K,  (1  -  0)| 


4  9(:. 

<sySS*exp(jas  -  ja*s  )>  *  -  j 2 w  a,  a  —  Mk,  -  k% )  6  (ky  -  ky ) 

•  W  (k^  +  k  sinfl,  ky  )  exp  tt  °\  la  +  a*2  -  2aa*0, 

I  9’(:. 

<Vy  SS*  exp(jas  -  ja*«’) >*  -  2tr6(kx  -  kx )  8  (ky  -  ky)o,  •  W(kx  +  k  smO.k^ ) 


*  A  -•©  (£)]  ■ «p[  i  A  i*’ + ■*’  ■  2"'c. '} 

[j*C 

7/ 

+ °? M#  (it)  ]  cxp  j-  i  °i  i*a  +  a#a  -  2aa#c,  i| . 

W(kx,ky)  is  the  roughness  spectrum  of  s(x,y)  and  is  equal  to  the  two-dimensional  Fourier 
transform  of  Ct  (u,v)  .  Usinf»  the  form  of  Cj  (u,v)  previously  given,  Wfk^  +  ksinfl,ky)  = 

<IhhIhh>  and  the  integration  is  carried  out  in  k  and  k  ,  the  delta  functions  will 
become  one.  Placing  the  expression  given  previously  for  C(  (u,v)  and  C(u,v)  in  the  aver¬ 
ages  and  changing  the  integral  to  polar  coordinates  results  in 


<WhV  =*K* 


=  KK*////{  2  no]  A  A*  +  (A#B  +  B*A)  8ir  jk  cost)  o 2 

+  j(AC*a  +  A*Ca*)o*  'i'"*  e' +  (A*D  +  AD*)  . 
jHk  cost?  a  r  cosa  n  of 


r  sina  a. 


+  (AE*a  +  A*Ea*)  j4ir 


0*  r  cosoe"r^  4rro,  o  BB*  B*  °2  BB*K,  r2  sin2  a 


-  l6rr(BC*a  +  B*Ca*)  k  cos 0  o1 


1  ■  j  4  -r2/lJ 
r  sin  ao.  e 


Biro2  o2(B*l)  +  BI)*)K  r2sinacosa 

— 1 - ■■4,,  a  -  -  1 6ir  0*  o  (BE*  a  +B*Ea*) 

I  .  (  a 


3,.3 

kcos0r  sino  cosoe 


Xital'l'*  „  j  .  3  3 «, 3 

+ - jr~  —  1  -  -2r  *jn--  e'r  /!  -  I6»r  (CI)*a*  +  C*l)a) 
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.  -ii.  4  -r  /I 

kcos®o  r  sinacosao,  e 


I2L2G 


(CE*  +  C#E)8wa4  r2sinacosae 

r5 


*'rJ/,a  4iro*  c2  DD* 

” +  lTg 


3/|l 

8»ro2  a2  DD#K,  r2  cos2  a  I6ir o*  r2  cos2  aka2  cos 0e  r 

1  1 - (DE#a  +  D*Ea*) - - 


W 


I'L'G 


+  - 


i,.i 

4rr  EE#o4|  e*r  7  8»rEE*a 


i41r2cos2ac*r,/|,\ 

l4  / 


W(k  +  k  sin0,k  ) . 


•  «xp{-  J  °\  (a*  +  a#,)|  I1  +  °t  aa*C,|exp  [-  K,  r^o]  * 


exp(-2jkr  sin®  sina)  .  rdky  dkydrda  . 

In  the  above  expression,  the  term  exp  (o2  aa*C( )  has  been  replaced  with  I  +  o^aa'Cj  . 
The  final  form  of  <IhhIhh>  must  be  left  as  a  double  integral  in  k^  and  ky  .  Integra¬ 
tions  can  be  performed  in  a  and  r  as  with  the  quasi- specular  term.  In  order  to  do  this, 
a  few  more  integrals  must  be  specified: 


/ 

o 

/ 


cosa  exp(-j2kr  sin®  sino)  da  =  0 


cosasina  exp(-j2kr  sin®  sina)  do  =  0 


m 

fCi 


a  ,  ...  ,  .  2w  J.  (2kr  sin®) 

cos  a  cxp(-i2kr  sin®  sina)  da  =  - '■ - 

H  J  2k  r  sin® 


The  integration  in  r  can  be  carried  out  using  the  same  Hrssel  function  integral  that  was 
used  earlier  to  evaluate  the  quasi-specular  term.  When  the  integration  is  carried  out  in 
•  and  r,  the  following  equation  results  for  : 


<1  I  *  >  = 
*HH  *HH 


m  am 

ii  ff  (  AA#L2Gexp|-k2L2sin2 

kk‘M  ( — ^ — 


«g/kj 
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I 

2(BC*a  +  B*Ca#)kco80o’  o’  rxp  |  -k’  sin’  oj  jp™  +  prj 


4(BC*a  +  B*Ca*)k’cosfl  sin’flo’  o’  fxp/-k’sin’ 


L’GlTA  *  4-1 

Li’g  i’J 


o’CC*L’gi’ 


■  2  .  1  „ 
-k  sin  0 


+  (K  I’  +  i’<;7  "xp  ,  » 

LvTTo  t 


*  |  2  2  _ 
K  +  2L  Cj  k  sin 

1  (K,  I*  +  L’G) 


2(l)K*a  +  l)*Ka*)o’  o’  k  «:os0  rxp  ^  -k’sin’  +  ’jT'J 

Hvtf 

o|aa*AA*rXp|-k’«in'y,|j^:  4  pjl 

2[i^  *  f] 


2(AB#  +  A#B)aa#a*  0*^0060  sin0  exp' 

,  i  .  1 . 

l  -  k  sin  0y 

L’g[^ 

Ll’g  i*. 

a 

+ 


aa*BB#aJo,L,GI* 
(K,  ?  +  L3  G)2 


**xp 


1 1  .  «. 

•  k  sin  6 

L*G  +  1* 

1  + 


2k,ginJ0K|  l4 
(K,  I*  +  L*G) 


+ 


aa'DD'ojVL’GI1 
(K,l*  +L,G)J 


Terms  containing  o*  have  been  dropped.  In  order  to  obtain  a  numerical  result  for  a"  H 
as  a  function  of  0 ,  the  double  integral  in  must  be  evaluated  numerically. 

b.  Depolarized  Term.  Attention  is  now  directed  to  the  derivation  of  the  depolar¬ 
ized  scattering  coefficient  o°  v  .  The  basic  expression  needed  to  compute  o“  v  is 


o  _  lim 

7hv  ft  - 


4*R’ <E„vE»*v> 


e,e; 


The  objective  now  is  to  calculate  the  quantity  <EHyEjv>  .  The  i  and  k  components 
from  equation  (2)  form  EHV  : 

Ehv  =  -  K  |l  +  Ey  Zy  (*sc 0  cos 0  +  s^  Ej  -  rjlljf  ?•„  csc0  cos30 
- 17 H—  Zy  sin0  -  r?llj  sy  sin0  +  tjII_  co t>0  +  17 ll_  L%  sin0 

+  rjll_  s^sinfl  +  nllj  sycosd|  exp  (-jk  n(  -r)dydx  . 

Placing  the  surface  fields  from  equations  (4)  into  the  above  expression  for  EHV  results 
in  the  following  equation: 
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•jZy  +bj6y|  exp  [- 2jkx  sin0  +  2jkZ  cos0] dydx  +  I, 

where:  a,  =  2R  cse0  cos 0;  b,  =  2R  sin0  cos0 
2  0  2  0 


HV 


■  ‘Mi  A'  +  B'Zx  +Csx  +  L)Zy  +  £’sy|  S  exp(jas)  exp  (-2jk  x  sin  0 


2jkZcos0  ]  dkx  dky  dy  dx 


A  =  -Q(l  +  R„)  ky  kx  -  Q(l  +  Rn)  cosfl  (k, -k,) 

B'  =  -  Qb ky  kx - Q g cosfl  (k,  -  k^ )  -  Q  ( I  +  R0 )  sinO  ip.(k,-k't) 

k  k 

C'  =  Q(1  +R0)kyk'x  -  Q(1  +  Ro)sin0  (kx  -k'x) 

[k  (1^  +  k  k' ) 
— - — 


r  k  (k1  +  k  k  >  ka  k’  i  r  k  (kx  +  k  k") +  kj  k 

g«(i  +  R,)0-i"«p-%-k  —  +  -Vj'9(1+R°)l  k 

It  should  be  noted  that  IHV  has  been  written  in  the  same  form  as  lHH  and  that  only  the 
definitions  of  the  coefficients  of  the  slope  terms  are  different: 

=  KK.//7/ Z, t«,b1Zy»;  +»,l):Z/y  ♦b'.y,} 

•  exp  j-  2jk  sinflfx  -  x  )  +  2jk  cos0(Z  -  Z’)|  dy  dxdy’dx’+  ®hv*jTv  • 
In  order  to  compute  <EHV  EH*V>  ,  a  few  more  averages  are  needed: 

<ZyZy  •  xp<2,k  ,  -  ■/.))>  --  -a’  K,  (£)]<■«,,{.  K,(l  -C)| 


CO80 
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<»r>  =  <s;>  =  o 

i  **c. 

.  it1  • 

Using  the  above  averages,  the  equation  for  C,  (u,v)  and  the  approximate  equation  for 
C(u,v),  the  expression  for  <  E^y  EH*V  >  becomes 


<EhvEh*v>  =  KK 


f  f  (2a  a?  4K,  o*  r*a?  cos* 
'{  \  L’G  L4Gj 

4 o|  r2hJoosJae  r  ^  1 
- p - j  exp(-2jkr  sinfl  sii 


a, 

„  i .  a  -r  /I 
a  2o.  b,  e 
-  +  — - 


exp(-2jkr  sinfl  sina)  exp[-K,  (1  -C)|  rdrda 


+  <IHV,H#V>  * 

The  change  from  rectangular  coordinates  (u,v)  to  polar  coordinates  (r,a)  has  been  made 
in  the  integral.  The  integrations  in  a  and  r  can  be  carried  out  easily  hy  use  of  the  pre¬ 
viously  defined  integrals  and  the  approximate  expression  for  C(r).  When  these  integra¬ 
tions  are  performed,  the  final  expression  for  <EHV  E*v>  is 


2irKK*hJ  a,2  Kj  l^GI1 

(Icjl^+T/Gy1 


^  ^HV  *HV  > 


The  equation  for  <‘Ihv*hv>  ^  the  same  as  that  for  except  that 

A,  B,  C,  I),  E  will  be  replaced  by  A,  B,  C,  I),  E,  respertively.  It  is  interesting  to 
see  that,  although  the  depolarized  scattering  coefficient  does  depend  on  the  large  sur¬ 
face  undulations,  if  o,  =  0  ,  the  scattering  coefficient  goes  to  zero  also.  A  discus¬ 
sion  of  the  results  and  their  meaning  will  tie  given  in  a  later  section  of  the  report. 


5.  Vertical  Polarization. 


a.  Like-Polarized  Term.  Consider  now  a  vertically  polarized  plane  wave  of 
unit  amplitude  with  a  time  harmonic  exp(jut)  incident  upon  the  composite  rough  sur¬ 
face  p(x,y)  =  Z(x,y)  +  s(x,y) .  The  surface  is  again  assumed  to  be  a  dielectric  with  a 
gaussian  distribution  of  surface  heights.  It  will  be  required  to  derive  expressions  for 
Oyy  ,  the  like-polarized  backscatter  coefficient,  and  o°H  ,  the  depolarized  backscatter 
coefficient.  In  order  to  do  this,  expressions  must  be  derived  for  Eyv  ,<Eyv  E*v>  , 
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EyH ,  and  <EyH  E^  >  ,  the  fields  and  mean  power  densities  of  the  like-polarized  and 
depolarized  returns.  Then,  o£y  and  can  be  determined  by  the  following: 


o  lim  4»R*<Evy  E^y> 
°W  R"li  £  £# 


o  _  lim  4wRj<EvhE*h> 
'Vh  “  R-«-  — 


E|E* 


For  a  vertically  polarized  wave,  the  incident  electric  field  vector  is  written  as 

*•  4  ^ 

E,  =  -  (icosfl  +  ksind)  exp  [-jk(xsind  -  zcos 6)  J  . 

From  Maxwell’s  equation,  the  incident  magnetic  field  is 


H  -  JCXP  {-  jk(xsinfl  -  z  cosfl)} 


where  r?  is  the  intrinsic  impedence  of  free  space.  For  vertical  polarization,  it  will  be 
easier  to  work  with  the  magnetic  field  vector  rather  than  the  electric  field  vector.  The 
locally  incident  magnetic  field  vector  H,'  can  be  written  in  terms  of  the  local  coordi¬ 
nates  (X,y,T)  in  the  same  manner  as  the  incident  electric  vector  was  written  in  terms 
of  local  coordinates  for  horizontal  polarization: 

H]  =  i(x<  x  •  j  )  +  y(y  •  j )  +z  (z  •  j)Jexp|- jk(x  sin#  -  Zcos0)j 
*  exp|-jk[n,  •  x ) x  +  (n,  •  y)y  +  (n, 

The  above  vector  dot  products  have  been  determined  previously  and  are  here  repeated 
for  convenience: 


x  •  j - Zy  D0  cos0 

-  i  _  , 
y  -j  =  1 

z  -  j  =  -  Zy  . 

Once  again,  since  J  •  y  =  1  the  effects  of  the  x  •  j  and  z  •  j  components  will  be 
neglected.  Using  the  quantities  n(  •  x,  n,  •  y,  and  ri(  -z  as  computed  previously,  the 
local  incident  magnetic  field  vector  can  be  written  as 
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H|  »  ^  y  exp  jk(x  und  -  Zcos0)|exp  j-jk(3Tsin0-Y  cos0  )| . 

The  problem  locally  near  the  surface  is  that  of  a  vertically  polarized  wave  incident  onto 
a  slightly  rough  dielectric  surface  with  an  angle  of  incidence  equal  to  $'  .  This  problem 
is  worked  in  Appendix  A,  and  the  results  there  are  restated  here  in  the  local  coordinate  sys¬ 
tem.  The  fields  are  correct  up  to  first-order  terms: 


f  ° •  *»  \ 

HI=  '  // B,1(fc,.k,)EYPd  k,dk,[ 

>  .»  .M  / 


EXP 


H-  =  | cxp(-jkx  sinfl  )  [exp(jkz  co80)  J  R||  exp(-jkTco80)] 


09  00 


ky )  EXP  dks  dky  j  EXP 


Hi  »  I  ff  ?,<k,.k,)EXPdk,dk(jEXP 
%  *t\ff  IMm >  * k1G„(k.,kr)|EXFdk 


-  exp(-jkx  8in0)|kco80'exp(jklco8d)  -  kRn  co80'exp(-jkTicos0)l  »  EXP 


} 


/  /  |k1l>„<k,,ky)tkIK„(kI,k,)|  EXPdk.dk  )  EXP 


,  / 

Ej  =  Y  \  sin0exp(-jkx  sin0  )  (exp(jkz  coe0  )  +  exp(-jkz  cos0  )) 

+  ff  Ik*Gy(k*’ky)-kvD*«(k«  kr))!:XPdk*dkyJ  EXP 


EXP  =  exp|  +jkK  x  +  jky  y  -  jkj  ] 
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EXP  =  —  exp  [-jk(x  sin0  -  Z  cos0)j 


d.A-V  - 


A„  a,,  -Ba.- 

GyA»ky>  =  i~Wi  a  1 

"ii  “za  "ia  *ai 


F  ,k  ■■  >  -  WW  +  WW 
*  u  *  *  *  S  *  I, 

*1 

Ao  =-i,k\ak#Q1  -  a?kfk2ky  Q,  -W^kk*  -  Vk^kk'5 
Bo  =kxk;kk'JSo  -  uk’k^k.Q, 


So  =  jky  S(kx  +  k  sin  0’,  ky  )r?T||  sin0'(kJ  -  k*  )/2wk'7 

V  =  - jT||i7kaina0'(k*  -k'a)S(kj|  +  ksin0,’ky)/2irk’a 

Wo  =  j(kx  +ksin0)ijT||  8in0’(ka  -k'2)S(kx  +  k  sin0,  ky  )/2wk 

Q,  =  j(kJ  -  k'2 )( 1  -  R|j)  cos0'S(k|(  +  k  sin0’,  ky)/2irk 
an  =  *7kykak'a  +rjkiJkxk  J  +  ijkJl^kf  +  TjkJk'akt 

*0  =  Tjkykik  *k*  +  ^  kyk*k* 

aJ(  =T?k'ikIk'tky  +  rjk,k>ikyki 


-  ijk’1  k’k',  +T?k’Jkakx  +  ijk'k^'J  +r?kJkxkl  , 

where  R|j  is  the  local  Fresnel  reflection  coefficient  for  a  vertically  polarized  wave.  The 
quantity  Tj|  U  equal  to  1  +  R^|.  The  form  of  R||  is  as  follows. 

i  rr*  .a  i v 


a.  . .  i  a.  a. 


k’  CO80*  -  k*/  k'  -  k  sin  0 


k  ^cos©’  +  k  /k’ 


Bine’  *  sinff  -  L%  cose 


cose'  •»  cose  +  Zx  sin  $  . 


In  order  to  obtain  a  workable  expression  for  P  ,, ,  it  will  be  necessary  to  approximate  Rj, 
by  the  first  two  terms  of  its  Taylor  series  expansion  in  Zx ,  about  Zj  =  0.  When  this 
is  done,  Ry  becomes 


k'cos  e  -  k^/k'1  -  k*  ain^'e 
^  k’*  cose  +  k  y/  k’1  -  k*  sin*  d 


t  =  2k'ak  sine  tk*a-  k*l 

^k^?s!n*T"{"k'1coee  +  k^k’a-  k*  sin^  } 

At  the  surface  where  the  fields  are  to  be  computed,  x  =  0,  y  =  0,  and  z  =  §(x,y) .  It  is 
also  necessary  to  write  Dx| ,  Gy| ,  and  Ffl  in  terms  of  slopes  so  as  to  make  a  manageable 
solution  later  on: 

DKi(k*’ky>  =  ici  <di  +  d2  Z,)S<k*  +  k«I'^ky) 

Gyl<k*’ky)  =  jCl  (8,  +KjZ*)S<kx  +k8in0'’ky) 

Fn  (k*’ky)  =  ici  (fi  +f2  z,)s<k,  +ksinfl',ky) 

c«  =  f“i.  aaa  "  a.2«2. 1"‘ 


dl  *11^2^0  "  *21  ®8 

dj  -  *||  "  *j| 

% 

8|  ~  a8*a  "^2^0*12 

^2  “  °9a2J  “^2^1  *B 

fi  =  <ky*»  +  Ml>/k. 
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fl  =  <ky*1+ksd3>/k, 

«,  =  _ro)(a3  +«t)  +  °5<1  ♦r0)(k,«n#  +  ksin,0)  +o4(l  +r0) 

0,  =  -«,rc  +0,  (1  -  r0)  -  a,r^  +  a4(l  -  r0)  -  a,  (1  +  rc)  <2ksinficosfi  +  k^cosff) 
+  <Vo*k«8in®  +k«n,fl)  +  atr^  +0,(1  +  r0) 

0O  =  k*(1  +  »«)  “n*  ~  k,<1  -  r0)cosb 

0,  =  k« r„anfl  -  kj(l  +  rc)  cosfl  -  ks(l  -  r„)  sin#  +  k^costf 


k  nkjkCk1-  k'a) 

**  si 

o,  =  -ijl^kki(kJ-k'J)co8J/2(r 

o2  =  -tjkk^kj(kI  -  k'a)sin0/2ir 

a3  -  -  rjk^k^k*  -  k'J)  cos0/2» 

a4  =  -rjkikk^(k1  -  k*)  sin0/2tr 

Qj  =  -nk^kfk1  -  k’2)/2 jf 

a6  =  t?kik'ikJ  8in,fl(k*  -  k’a)/2ir 

a,  -■  -  2tjkik’l  kJ  sinfl  cos0(k*  -k’a)/2r  . 

In  deriving  the  above  expressions,  terms  containing  Za ,  Z* . etc.,  have  been  neglected. 

The  electric  and  magnetic  fields  evaluated  on  the  surface  become: 


D„«P<-jMdk*dky 


} 


EXP 
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H-  =  jexp(jk»  cos®)  +  (r#  +  rjZx  )  exp(-jks  cos®)  +  //.  yl«p(-jkt8) 
dk„  dky|  EXP 

Hi  =|  ffYn  exP(ikI*)dk,  dky|  EXP 
h  =  ff  <kyF,.  +  k.Gy.>exP<*ik,8>dk«dky 

\.M  .M 

-kZxsin®  [exp(jksco80)  -  (rD  +  r0'Zx)  exp  (- jks  cos®)) 

-  k  cos®  exp(jks  cos®) +  k  cos®  (r0  +  roZx)exp(- jks  cos®)\  EXP 


+  kx  Ffl )  exp<-jk 


..Hk.dk,  J 


Ej  =  -jj-  j-k(sin®  -  Zxcos0)[exp(jks  cos®)  +  (rc  +  r^Z,  )exp(-jks  cos®)) 

•//.  kG»I  ~  ky  D«l  >  CXP  Jk«S>  dk*  dky  }  KXP  • 

-  mm  .m  ' 


(5) 


The  (i  cos®  +  ksin®)  term  in  equation  (2)  will  now  form  the  like-polarized  backscattered 
field  when  the  incident  wave  is  vertically  polarized.  Using  this  term  from  equation  (2)  and 
the  above-defined  surface  fields,  we  can  obtain  an  expression  for  the  like-polarized  back- 


scattered  field  Evv : 


Evv  =  -  K  f  f  ( a’(  +  bj  Zx )  exp|-  2jkx  sin®  +  2jkz  cos®|  dydx  +  lvv 


a’(  =  2rocos0  bj  =  2rosin®  +  2r„cos  0 
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•w  -  K /// /  {\  ♦  B,  Z,  +  C,  \  +  D,  Z,  ♦  E,  v  }  S  «PQ..) 

•  exp  |-2jkx  sin 6  +  2jkZ  cosfl  jdkx  dky  dydx 
S  =  S(k8  +k«ind',k, )  a  =  kcosfl  -  kf 


\  =  {  kr  +  k«^i  +  k  co>®} 

B.  =-j“  {kyf2  +  ktgj  +  gjkcosff  +  g,k  sinj| 

C.  =-^r{k«*.  ‘  krd,+8ik  ™9  } 

D.  =rab  {  (M. +k.f,)cOB«  +  kd.} 


f,  cosfl  -  dj  sin0 


Use  has  been  nude  of  the  derived  forms  of  D^.C  and  .  It  can  easily  be  seen 
that  Eyy  consists  of  two  terms  the  first  of  which  depends  solely  on  the  large  undula¬ 
tions  and  represents  the  quasi-specular  component.  The  second  term  Iyv  ,  although  it 
depends  on  the  large  undulations  through  two  slope  terms  and  a  phase  term,  arises  be¬ 
cause  of  the  slightly  rough  surface  s(x,y).  The  form  of  Evv  is  now  the  same  as  that  for 
Ehh  ,  the  only  difference  being  in  the  definitions  of  the  coefficients  on  the  slope  terms. 
An  expression  for  <EyV  01,1  ^,Cn  wr'ttcn  directly  without  further  work: 


<Ew  E£y> 


=  KK* 


G  8s a  L’Gs',  I)’  kJ  sinO  cos0 


K 


4sL,Gb|1’,<7J  k* 

KJ 


—  —  ^  exp  |-kJ  L* G  sin^/K^ 


+  <IwIvv> 


The  equation  for  <Iyy  I^v>  will  be  the  same  as  that  for  <IHH  I  *H>  except  that  A, 

B,  C,  D,  and  E  will  be  replaced  by  A, ,  B, ,  C, ,  D( ,  and  E,  respectively.  Once  again,  it 
can  be  seen  that  <lyvI*v>  will  disappear  when  o(  =  0  leaving  only  the  quasi-specular 
term. 
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b.  Depolarised  Term.  The  depolarised  ecettering  coefficient  for  a  vertically 
polarised  incident  wave  is  <fVH  .  The  objective  now  is  to  obtain  an  expression  for 
<Evh  E*|>  .  The  J  component  of  equation  (2)  will  give  s  basic  equation  for  EyH 
in  terms  of  the  surface  fields: 


Evh=  * 


E-  Zy  esc®  cos*®  -  Ey  cos®  -  Ej  ^  cos®  +  tjH- 


+  ijH^  Zycsc®  cos®  +  ^  +  EjZy  sin®  +  ^  Ej  sin® 

-  E^Zasin®  -  sa  Ep  sin®  jexpf-jkn,  *r)dydx 

The  pr  eviously  stated  surface  fields  used  in  the  above  expression  for  EVH  results  in 
the  following  equation: 


+  bjs  l  exp[-2jkx  sin®  +  2jkZ  cos®)  dydx  +  ly||  , 


where 


aj  =  2r0csc®  cos®  bj  =  2rc  sin®  cos® 

I*-*////  *  Bi Z.  +  CA  *  Di  zr  *  Ei  \}  S  «P(j*.) 

•  exp  (-  2jkx  sin  0  +2jkZ  cos®  ]  dks  dky  dydx 

(kid|cos®  +  ks  fj  cos®  +  dt  k] 

B',  =  X  KMi  +  k,  fj )  co*9  +  d,  h  +  (kidJ  +  kE  f, )  sin  0 ) 

C*  =1T  {(k,di  +kEf1)sinfl+f|k| 

Di =  nb  lkvf.  +Mi +k*.  COtt^ 

Ei  =-Y  l(k,f|  +klg1)sinfl  -(kEgj  -  kyd,  )cos®J  . 


I 


f 

1 
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The  depolarized  expression  consists  of  the  aim  of  two  terms,  the  ftrat  of  which  ia  depen* 
dent  only  upon  die  two  y  slope  terma.  The  term  IVH  ia  similar  to  expreaaiona  derived 
previously.  The  form  of  EyH  can  be  seen  to  be  the  same  as  that  for  E^y  so  that  an 
equation  for  <EyHE^H>  can  be  written  down  directly: 


<EVH*VH> 


2tr  KK#bj*  oJK,  L*  Gl* 
(K,l*  +  L*G)3 


exp 


,1.1. 
-k  am  0 


+-T 

lg  r 


+  <IVHIVH> 


The  equation  for  <IVH  I  *H>  isof  the  same  form  as  <IHHI*H>  except  that  A,  B,  C, 
D,  and  E  will  be  replaced  with  A', ,  Bj ,  Cj ,  Dj  and  Ej  .  Again,  the  depolarized  term  dis¬ 
appears  if  o,  =  0 ,  although  the  expression  is  still  a  function  of  the  large  undulation 
parameters. 


6.  Consideration  of  Local  Vertical  Components.  ?•  hen  the  incident  fields  were 
written  in  the  form  of  local  coordinate  components,  the  effects  of  the  x  *  j  and  z  •  j 
components  were  neglected.  It  is  the  purpose  of  this  section  of  the  report  to  calculate 
the  changes  that  would  occur  in  the  resultant  backsrutter  coefficient  if  the  local  vertical 
terms  are  not  ignored.  The  changes  needed  for  an  incident  wave  with  horizontal  polari¬ 
zation  will  be  worked  first.  For  a  horizontally  polarized  incident  wave,  the  incident  field 
was  given  earlier  and  is  repeated  here  for  convenience: 

E,1  =  |x(x-j)  +  y(y  -J)  +  z(z  -J)j  exp  j- jk(xrin®'-z  cos®')  j  EXP 

EXP  =  exp  [-  jk(xsin®  -  Zcos®)] 

* 

x  .  j  =  -  Zy  D0  cos  0 

y  j  -  i 

z  •  y  =  -  Zy 

Do  =  (sin®  -  Zxcos®)” 1  . 

The  amplitude  of  the  local  vertical  component  can  be  computed  from  Maxwell's  equa¬ 
tion  written  in  local  coordinates: 


jVxE? 

k 


The  problem  is  non  to  compute  the  amplitude  of  the  y  component  of  H|  .  This  am¬ 
plitude  will  then  be  multiplied  by  each  of  the  surface  field  components  determined  for 
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vertical  polarization  and  then  added  to  the  aurface  fields  already  given  by  equations  (4) 
for  horizontal  polarization  in  order  to  obtain  the  new  expressions  for  die  total  field. 
The  y  component  of  H|  will  be  written  as  H|-  and  is  computed  to  be 

Z  esc  6  /  \ 

Hj-  *  — ^ —  exp  <-jk(xsin0-Tcoa0)  > 

In  this  calculation  the  higher  order  slope  terms  have  been  ignored.  The  total  fields  on 
the  surface  can  now  be  written  in  the  form: 


Ej  -  E|T  +  Ejj 


E-  =  E  -  +  E„- 

^  ly  iy 


Ej  =  ElT  +  EJ? 


Hi  =  H.I  +  Hai 


Hy  =  H,7  +  Hjy 


Hr  = 


The  subscript  1  refers  to  the  surface  field  determined  by  a  local  incident  horizontally 
polarized  wave,  while  the  subscript  2  refers  to  the  surface  fields  determined  by  a  local 
incident  vertically  polarized  wave.  The  surface  fields  due  to  a  local  incident  vertically 
polarized  wave  can  be  determined  by  multiplying  Zy  csc0/tj  by  the  fields  calculated  in 
Appendix  A  and  evaluating  on  the  surface: 

mm  m 

Z  C8C0  (  f  f 

Ei*  =  — I —yJJ  <k,F«.  +  k,G,.>e*PHMdMky 


E* 


-  kZx  sin0  [exp  (jks  cos0)  -  (rc  +  rn  Z% )  exp  (- jks  cos0)] 

-  kcos0  exp(jkscos0)  +  kcos0(ro  +  r^Z^)  exp  (-jkscos0)^  EXP 

-Z„ 


■} 


=  (k«D**  +k*F*»)exP('jk.8)dk«dkyJ 


EXP 


Z„  esc  6 
=  y 


_ - - yy  j  (k„GyI -kyI)xI)exp(-jkls)dkBdkjr 

-  k(sin0  -  Zi  cos0)  [exp  (jks  cos0)  +  (rc  +  r„  Z^)  exp  (-jks  cos0) )  VEXP 


tjHj;  =Zycsc 


{//» 


,i  exP(‘jkgs)dkiidky  p  EXP 
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■> 


r?Ha-  =  Zy esc®  | exp(jks  cos®)  +  (r0  +  r*  Z, )  exp  (- jks  cob®) 

•»  m 

*  f  f  G(1«p(-jk1.)dk,dky}EXP 


exp  (-jkfB)  dksdky 


} 


EXP  . 


The  above  fields  must  be  added  to  the  fields  given  by  equations  (4)  in  order  to  obtain 
the  total  surface  fields.  The  above  equations  are  correction  terms  owing  to  the  fact  that 
local  vertical  polarization  components  exist.  When  the  total  surface  fields  are  now  sub¬ 
stituted  into  the  expression  for  Em ,  it  is  found  that  the  quasi-specular  term  remains 
unchanged  if  high-order  slope  terms  are  still  ignored.  The  only  change  in  IHH  comes  in 
a  new  definition  of  the  D  term.  It  can  be  shown  that  D  now  becomes 


D  =  esc®  cos*®  kykxQ(l  +  R,,)  +  co8gkc8C?  (kid,  +  k^f, )c,j 

k  k 

+  jct  dj  CSC®  +  esc®  cos®(l  +  R,,)  — (kf  -  k' )Q  +  sin®(l  +  R„)kykEQ  . 

Two  more  terms  have  been  added  to  D  .  When  the  new  total  surface  fields  are  placed 
in  the  expression  for  EHV  ,  the  definition  of  a3  changes  and  the  definition  of  D’ 
changes  as  follows: 


D'  =  -  ^  (kyf,  +  kfg,  )jc,  -  esc®  cos®  jc, g3  +  esc®  cos®(l  +  R^MkjJ  +  kfkf)Q 


+  '0  +  R„) 


Mk’.+kjO  k’yk'f 
k  k 


Q  esc® 


a2  =  2  esc®  cos®  (r0  +  R„)  • 

The  term  a2  is  not  used  in  the  final  results  and,  therefore,  is  not  needed.  The  new  defi¬ 
nition  of  the  D'  term  has  two  more  terms  added  which  come  from  local  vertical  polari¬ 
zation  terms. 

A  similar  analysis  to  the  above  can  be  performed  for  the  case  of  vertical  polariza¬ 
tion.  This  time,  the  interest  lies  in  determining  the  amplitude  of  the  local  horizontal 
component.  This  can  be  done  by  using  the  following  Maxwell  equation: 
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v « A;  =jU«.  £; . 
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In  the  above  equation,  Hj  it  the  incident  magnetic  field  written  in  local  coordinates  for 
a  vertically  polarized  incident  wave.  The  equation  for  the  y  component  of  E| ,  E|-  ,  is 
then  1 

Eg  =  -Zy  csc0  exp  jk  (x  &;'d'  - 1  cos0)  j  EXP 

EXP  =  exp  j-  jk(xsin0  -  Zcos0)|  . 

The  equations  for  the  surface  fields  for  horizontal  polarization  must  be  multiplied  by 
the  amplitude  -  Zy  csc0  and  added  to  the  surface  fields  given  by  equations  (5)  to  ob¬ 
tain  the  total  surface  fields.  When  this  is  done,  the  total  surface  fields  are  then  placed 
in  the  expressions  for  Evv  and  EVH  to  determine  differences.  With  these  changes, 
the  definition  of  the  D,  term  which  is  the  coefficient  of  Zy  in  lyv  becomes 


D,  =jc,  |  —  8?j-c8cg  (dJkl  +  f,kx)  +  d,  8in9  +  d,cosJfl  csc0  j 

k  k 

+  csc0(l +  R0)kyk|Q  i  cscfl  cos*  (l +  RJ (kf-k’g)Q  . 

The  term  D,  is  the  only  change  that  occurs  in  Eyv .  The  definition  of  two  terms  changes 
in  EyH  .  The  terms  are  aj  and  Dj : 

a^  =  2  esc 0  cos0  (re  +  R^) 

jc  j 

Dj  =  -r-  [ky  fj  sinfl  +  kfg,  sin0  +  g,  k  csc0  cos0  +  kyf,  csc0  cos  0 

+  kgg(  csc0  co8*0  j  -  csc0  cos0(l  +  R0)(k*  +klk’f)Q 


-  esc 0  (1  +  R^) 


Mk2„  +  UO 


Q 


The  only  effective  difference  in  E^  is  in  the  D’(  term  since  a^  does  not  appear  in 
the  final  result. 


This  ends  the  analysis  portion  of  the  report,  the  next  section  discusses  the  results 
in  terms  of  graphs  of  the  backscatter  coefficient  versus  incidence  angle. 
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DL  DISCUSSION 


7.  Evaluation  of  Remits.  The  purpose  of  this  section  of  the  report  is  to  describe 
numerical  calculations  and  to  show  and  discuss  the  resultant  graphs  of  scattering  coeffi¬ 
cient  versus  incidence  angle.  These  graphs  represent  the  final  result  of  all  derivations  and 
are  plotted  for  different  surfaces,  different  soils,  different  moisture  contents,  and  two 
different  frequencies.  The  equations  for  ^HH ♦  o£v  ’  °vv  ’  and  °vh  involve  double  inte¬ 
grals  in  1^  and  ky  which  cannot  be  solved  easily  by  analytical  techniques.  It  is  neces¬ 
sary,  therefore,  to  use  a  numerical  integration  routine  in  order  to  obtain  final  calcula¬ 
tions.  A  two-dimensional  Simpson's  rule  wu  used  to  evaluate  the  double  integrals.  The 
question  of  limits  of  integration  for  the  computer  was  determined  by  using  the  exponent 
term  in  W(kB  +  k  sin0,  k  )  .  This  exponent  will  decrease  rapidly  as  k(  and  ky  in¬ 
crease.  If  e~ 10  is  assumed  to  be  insignificant,  then  the  limits  of  integration  (kw ,  kym ) 
will  be 

k..  =  k„  =  <2/5^ 

with  0  =  (f  .  A  computer  program  was  written  in  Fortran  IV  for  the  calculation  of  the 
backscatter  coefficients.  The  program  inputs  are  frequency,  limits  and  number  of  in¬ 
tervals  for  integration,  standard  deviation  of  the  large  and  small  undulations,  the  corre¬ 
lation  parameters,  and  dielectric  constant  of  the  surface.  The  output  of  the  program  is 
the  four  scattering  coefficients  in  decibels  for  different  angles  of  incidence.  The  func¬ 
tions  needed  for  integration  are  formed  in  four  separate  function  subprograms,  while 
integration  itself  is  performed  in  the  main  program.  Careful  examination  of  the  equa¬ 
tions  for  the  scattering  coefficients  shows  that  they  are  not  a  function  of  the  absolute 
values  of  a  or  L  but  only  their  ratio.  For  all  calculations,  the  value  of  a  read  into 
the  program  was  left  at  1.0  and  the  value  of  L  was  varied.  This  is  not  the  case,  how¬ 
ever,  with  the  small  surface  perturbations  as  both  o,  and  I  appear  to  have  importance 
by  themselves  and  so  were  varied  individually  in  the  program.  It  is  not  possible  to  eval¬ 
uate  the  scattering  coefficients  at  0  =  0°  due  to  the  csc0  in  the  1)1)*  terms  of  the 
double  integral.  The  cscd  term  came  about  because  of  the  assumptions  that  were  made 
in  the  derivations  and  these  assumptions  are  invalid  at  0  =  0°.  Computations  were 
started  at  0  =  10°  and  were  performed  in  10°  increments  up  to  0  =  80°.  Cal¬ 
culations  were  also  made  for  the  case  where  o(  =  0  or  when  the  surface  is  just 
a  smoothly  undulating  one.  In  this  case,  only  the  like-polarized  scattering  coeffi¬ 
cients  have  values  as  average  power  for  the  depolarized  terms  disappear.  These 
graphs  allow  a  determination  to  be  made  of  the  effect  of  the  small  surface  per¬ 
turbations  on  the  resultant  backscatter.  Two  different  frequencies  were  used  in 
the  calculations:  9.375  GHz  and  5.87  GHz.  For  each  separate  frequency  used,  it 
was  important  to  determine  the  validity  of  the  very  rough  surface  condition  by  evaluat¬ 
ing  K,  =  4kV  cos30  and  making  sure  that  K,  »  i  .  If  o  is  held  constant  (c  =  1) 
and  K,  is  evaluated  at  0  =  80°,  then  only  the  lower  frequency  need  be  calculated  as¬ 
suming  it  satisfies  the  condition  K(  »  1: 
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K,  =  1823  when  f  *  5.87  GHz,  0  -  80° 

It  can  be  teen  that  the  very  rough  surface  condition  is  good  for  all  calculations.  When 
the  frequency  was  changed,  the  parameters  o,  and  1  were  also  changed  so  that  ko 
and  kl  remained  constant.  The  dielectric  constants  used  were  taken  from  Lundien1' 
and  represent  soils  with  certain  moisture  contents  and  dry  densities.  The  soil  type, 
moisture  content,  and  dry  density  of  the  soil  are  written  on  the  curves  along  with  the 
other  data.  The  graphs  which  follow  this  section  are  labeled  I  through  30  with  all  ap¬ 
propriate  data  labeled  on  each  one.  Preceding  the  o°  curves  is  a  curve  of  the  correla¬ 
tion  function  p(r).  It  can  be  seen  that  the  correlation  distance  is  larger  than  L.  By  far, 
the  most  important  factor  affecting  the  results,  at  least  for  like-polarized  terms,  is  the 
ratio  o/L.  The  effects  of  polarization  are  noticed  more  at  large  angles  of  incidence  than 
at  small  angles.  This  agrees  very  well  with  the  experimental  results  discussed  by  Beck¬ 
mann  and  Spizzichino.12  This  would  indicate  that  if  surface  discrimination  is  wanted  by 
polarization  variation,  then  it  possibly  should  be  done  at  large  angles  of  incidence. 

The  cause  of  depolarization  is  the  slightly  rough  surface.  Yet,  if  a  slightly  rough 
surface  is  taken  by  itself  so  that  o  =  0;  then  no  depolarization  results  for  backscatter- 
at  least,  for  first-order  perturbation  terms.  Also,  the  type  of  depolarization  that  results 
from  the  large  undulations  by  themselves  (o(  =  0)  depends  on  the  sum  of  the  two 
Fresnel  reflection  coefficients  and  does  not  appear  in  the  final  results.  This  is  very  in¬ 
teresting  because  the  composite  surface  ends  up  with  depolarized  terms  due  to  the  fact 
that  the  local  surface  fields  determined  by  the  small  perturbation  technique  are  com¬ 
puted  in  a  local  coordinate  system  which  is  tilted  with  respect  to  the  original  coordi¬ 
nate  system.  This  means  that  the  depolarization  results  from  a  tilted,  slightly  rough 
surface.  It  can  be  seen  from  the  graphs  that  the  two  depolarized  terms  are,  in  general, 
not  equal  over  all  angles.  The  cause  of  this  could  be  either  the  result  of  the  approxima¬ 
tions  made  in  the  solution  or  there  really  is  a  difference  in  the  two  terms.  It  is  not 
possible  at  present  to  determine  which  case  is  correct. 

The  effects  of  each  individual  parameter  will  be  discussed  briefly.  The  effect 
of  differing  moisture  contents  can  be  seen  by  examining  curves  i  through  4  and  curves 
5  through  8.  The  higher  moisture  content  results  in  a  higher  dielectric  constant  which 
yields  a  little  higher  value  of  a°  for  each  angle  of  incidence  and  for  all  four  polariza¬ 
tions.  Curves  9  through  12  are  of  a  different  soil  type,  but  care  must  be  taken  in  com¬ 
paring  these  graphs  with  curves  1  through  4  because  the  moisture  contents  differ.  If 
the  difference  in  moisture  contents  could  be  taken  out  there  probably  would  be  very 
little  difference  in  the  two  sets  of  graphs.  The  effect  of  changing  the  o/L  ratio  can  be 


11 J.  R.  Lundien,  "Terrain  Analysis  by  Electromagnetic  Means,”  Technical  Report  No.  3-693,  Waterways  Experiment 
Station,  Vicksburg,  Mist. 

11  Beckmann  and  Spizzichino. 
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seen  by  comparing  curves  1  through  4  and  curves  13  through  16.  This  effect  is  particu¬ 
larly  noticeable  at  intermediate  angles  for  the  like-polarized  terms  where  the  difference 
becomes  very  dramatic.  The  ratio  o/L  can  be  seen  to  affect  the  like-polarized  terms 
much  more  than  the  depolarized  terms. 

The  influence  of  the  correlation  distance  (I)  of  the  slightly  rough  surface  can 
be  seen  by  comparing  curves  1  through  4  with  curves  17  through  20.  Increasing  the 
correlation  distance  I  has  the  result  of  bringing  the  o°  curves  down  at  large  angles  of 
incidence.  For  the  like-polarized  terms,  this  effect  becomes  noticeable  for  0  >  40°. 

In  the  depolarized  case,  the  whole  curve  is  lowered;  but  the  differences  still  become 
greater  at  larger  incidence  angles.  The  influence  of  oj  ,  the  standard  deviation  of  the 
slightly  rough  surface,  can  be  seen  by  comparing  curves  I  through  4  with  curves  21 
through  24,  and  curves  25  and  26.  This  parameter  has  little  effect  on  the  like-polarized 
term  until  angles  of  incidence  0  >  40°  are  reached.  The  effect  of  c,  on  the  depolar¬ 
ized  terms  covers  the  entire  range  of  incidence  angles,  and  the  difference  increases  a 
little  at  large  angles.  All  of  the  above  comparisons  have  been  made  at  the  frequency 
of  9.375  GHz.  It  would  not  be  reasonable  to  expect  that  exactly  the  same  effects  would 
occur  at  other  frequencies.  An  idea  of  frequency  differences  can  be  obtained  by  com¬ 
paring  curves  I  through  4  with  curves  27  through  30.  For  all  four  polarizations,  the 
X-band  frequency  gives  a  slightly  higher  value  of  o°  for  all  angles  of  incidence. 


IV.  CONCLUSIONS 

8.  General.  The  following  general  conclusions  can  be  made  from  the  analysis  and 
the  graphs. 


a.  A  correct  theory  for  explaining  radar  backscattcr  must  c  insider  a  compos¬ 
ite  surface  as  has  been  demonstrated  in  this  paper. 

b.  Depolarization  results  from  the  effects  of  the  tilted,  slightly  rough  surface. 

9.  Specific.  The  following  specific  conclusions  can  be  made  from  the  plotted  data 
at  the  frequency  9.375  GHz. 

a.  The  ratio  o/L  affects  the  like- polarized  terms  much  more  than  the  de¬ 
polarized  terms. 

b.  An  increase  in  the  moisture  content  of  a  soil  increases  o°  slightly  for  all 
angles  of  incidence  but  does  not  significantly  change  the  shape  of  o°  versus  6  curve. 
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c.  An  increase  in  the  correlation  distance  for  the  slightly  rough  surface  re¬ 
sults  in  e®  being  less  at  high  angles  of  incidence  (f  >  40*)  for  like-polarized  term.  For 
the  depolarized  terms,  o°  is  less  for  all  angles  of  incidence. 

d.  An  increase  in  o,  ,  the  standard  deviation  of  the  slightly  rough  surface  un¬ 
dulations,  results  in  a  higher  o°  .  In  the  c aae  of  the  like-polarized  terms,  a  higher  o°  be¬ 
comes  noticeable  for  9  >  40*.  In  the  case  of  the  depolarized  term,  a  higher  o°  is  no¬ 
ticeable  for  all  angles  of  incidence. 

e.  It  was  shown  from  the  graphs  that  the  changes  that  occur  in  o°  due  to  a 
variation  of  parameter  inputs  is  not  large  in  many  cases.  This  would  indicate  that  if  a 
radar  is  used  for  determining  o°  for  purposes  of  identification  and  discrimination  then 
a  very  accurate  calibration  would  be  necessary. 
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APPENDIX  A 


DERIVATION  OF  THE  SCATTERED  FIELDS  FROM  A  SLIGHTLY  ROUGH 
SURFACE  FOR  VERTICAL  POLARIZATION 

Consider  a  vertically  polarized  plane  wave  incident  onto  a  slightly  rough  surface  s(x,y) 
which  is  the  boundary  between  two  homogeneous  dielectric  mediums.  The  geometry  of 
the  problem  is  given  below.  The  surface  function  s(x,y)  represents  the  z  distances 
from  the  xy-plane  to  the  surface. 


Z 


A  right-handed  rectangular  Cartesian  coordinate  system  is  set  up  with  the  xy-plane  form¬ 
ing  the  average  value  of  s(x,y) .  The  xz-plane  coincides  with  the  plane  of  incidence.  The 
angle  0  is  the  angle  of  incidence.  The  conditions  that  must  be  made  upon  s(x,y)  are 
that  it  be  Fourier  transformable  and  that 


|ks(x,y)|  <  1,  where  k  =  2 ff/X  . 


The  propagation  constant  in  the  medium  z  <  s(x,y)  is  designated  k  and  is  equal  to 
o)yfJTf  .  The  incident  magnetic  field  is  in  the  +y  direction,  and  the  components  of  the 
total  magnetic  field  in  the  space  z  >  s(x,y)  can  be  written  as: 


H 


M  OO 

=  //U*(k«’k> 


)EXP, 


dk  dk 
*  » 


(la) 
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Hy  =  exp(- jkx  sin0)  [exp(jkz  cos0)  +  ^  exp(- jkz  cos0)] 


(lb) 


(lc) 


where  EXP(  =  exp(jkxx+jk  y-jkfz)  and  Rj|  is  the  Fresnel  reflection  coefficient  for 
a  vertically  polarized  wave,  'fhc  parameters  kx  and  ky  are  the  Fourier  transform  vari¬ 
ables  and  kx  is: 

kx  =  ^kJ  -  k2  -  k^  when  k2  >  +  k^ 

kf  =  -  j  ^/k2  +  k^  -  k2  when  k2  +  k‘y  >  k2 

The  total  magnetic  field  in  the  medium  z  <  s(x,y)  can  be  written  as 


K  =  ff  y)EXP,  dkxdky  (2a) 

Hy  =  T,|exp(-jk’x  sin0)exp(jk'zcos0)+yy* Gy  (kx,ky)  EXP2  dkxdky  (2b) 


(k*»ky)  EXPj  dkxdky 


(2c) 


where  EXP2  =  exp(jkxx  +  jkyy  +  jk^z)  and  T||  is  equal  to  I  +  Kj|  .  The  angle  0  is 
related  to  0  through  Snell’s  law.  The  parameter  k’x  is  defined  as  follows: 

k’  =  /k 2  -  k2  -  k2  ,  when  k 2  >  kx  +  k2 

Z  Y  "7  *7 

k’  =  -  j  /k2*"  +  k2  -  k’2  ,  when  k’2  <  k2  +  k2 

If  the  surface  were  flat  so  that  s(x,y)  =  0,  then  the  fields  that  would  Ik*  obtained  are 
those  of  equations  (I)  and  (2)  with  the  doublcd-integral  terms  eliminated.  The  tech¬ 
nique  of  small  perturbation  requires  that  the  amplitudes  Ik*  written  in  a  perturbation 


aerie*,  e.g.  Dx  =  DxJ  +  Dsj  ♦  Du  + . (where  the  dependence  upon  kx  and 

ky  has  been  omitted).  The  fields  near  the  interface  can  be  written  as  follows: 

=//{D- +<D-'ik',D-',t(Uu'jk'zD“  + . } 


exp(jkxx+jk  y)dkxdk 


(3a) 


«.=ff  (D«  -'H^l  ♦  j 


•  exP(jk„  x  +  jky  y)  dkx  dky  (3b) 

[| 2  2  2 
1  +jk(cos0 - 1  2°8 —  +  R||  -jR^kzcosfl 

^kVco^j  +  jj  L1+(Crl.jk,«C,1)  +  (c„.jk,,G„ 

•  M  .«  \ 

4 . } 


kVc 


exp(jkxx+jky  y)  dkxdk 


(3c) 


OO  Oo 

K  -  T||C\p(-jkx  8iri0)  [l  +jkzco80-—  +  J  J  |(^  +  (GyJ 

py3  +  ik,zGy2 — l~Tl  +  •  |  •exP(lk,  x+jkyy)dkx  dl^  (3d) 

OO  oo 

"■  11  *(F“  ■|k'zK")+(K‘>  ♦ . j 


cxp(jkxx  +  jkyy)dkxdk 


(3c) 


OO  OO 

=//  {k,  +  +  )k.<,  >  *  (f»  ♦  jk>K„  -  ‘v^)  ( . j 


exp(jkxx  +  jkyy)dkxdky 


(3f) 
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where  the  factors  exp(-jk(z)  and  exp  (jk'tz)  have  been  expanded  in  a  series  for  the  case 
where  z  js  small.  At  the  surface  where  z  =  s(x,y),  the  boundary  conditions  nx(H-H) 
and  n  x  (E  -  E')  can  be  written  in  component  form  as 


<H.  -  H>  (if)  <H.  -  h.> = 0  <«»> 

<»*-»;>  =°  •  (-“>) 


r,  f3H* 

bin 

1 

3H* 

3  Hy  1 

j 

W  VW  " 

dz  J 

H 

T7  1 

3y 

"37]  + 

I 

z 

=  ?(x,y) 

z  =  s(x,y) 

3Hx 

~w) 

k 

(>H, 
iTx  ■ 

■3)1 

=  0 

z  =  s(x,y) 

(4c) 

1 

jl[3H.  !!i 

jk  i  dz  “ 

L 

+  is] 

s(x,y)  k‘ 

p>r. 

-  1  dz 

sh; 
-  IT 

(4d) 

($[*(£ 

311 

*  "ay 

(M, 

\3» 

‘  “5y7 

- 

z  =  s(x,y) 

In  order  to  determine  the  six  unknown  amplitudes,  six  independent  equations  are  needed; 
so,  the  four  equations  above  and  the  two  divergence  conditions  V  *11=0  and  V  •  H  =  0 
combine  to  yield  the  necessary  number  of  equations.  The  partial  derivatives  can  be  ob¬ 
tained  from  equations  (3) .  When  this  is  done  and  only  the  first  order  terms  are  kept,  the 
following  results  are  obtained: 


k*  D*l  eXP(jkx  X  +  iky  y)  dk*  dky 


3z 


exp(jkxx  +  jkyy)  dkx  dky 


exp(jkxx  +jkyy)  dkx  dky 


(5a) 


j 

(Sl») 


(5c) 
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oo  a* 

7^=j  ff  Mi,  exp(jk,x+jkyy)dk4dky 

j 

j 

.3 

i 

<5d) 

M  M 

^ =  )  f f  kyGyi  «p(jkxx+jkyy)dkx  dky 

.  oo  .ae 

<5e) 

oo  oo 

^T  =  -j  / /  k;F;i«P(jkxx+jkyy)dkxdky 

•  oo  .  oo 

| 

(50 

1 

oo  oo 

Ijf  =  iff  kyFI.exp(jkxx+jkyy)dkxdky 

•  OO  .  oo 

I 

I 

(5g) 

OO  ny 

^ =iff  kxFi.  exP(jk„x+jkyy)dkxdky 

.00  .00 

(5K) 

dH 

=exp(-jkxsin0)|jkcoHfl(l-H||)-k2zco8,fl(l  +  Rj()j 

oo  00 

-iff  kiGy.exP(jk,x+jkyy)dkxdky 

-00  .OO 

(5i) 

3,,y  l  \ 

-jj-  =  -  jk  sinJ  cxp(-jkxHinO)  <  T(|  +  jkzco«0(l  -  R(|)  V 

OO  OO 

+  j  ff  kxGy|  '‘xp(jkx  X  +  jky y)dkx  dky 

* 00  .00 

<5j) 

on  00 

^7=-j  Jf^^W, 

-  00  .  00 

(5k) 

I 

00  00 

an  r  r 

3x  =  ij  J  KKie*p<iK*+fryy)dKdK 

-  00  .00 

(51) 

HO 


M;,«Ptik,x*jk,y)dk,dkr 

(5m) 

A  H  * 

=  T|j  exp(-jkx  sin0)  jjk  cos*  -  zk  cos  *  j 

+  j  jf  k;Gyl  exp(jky  x  +  jky  y)  dky  dky 

(5n) 

am  am 

£-j//  k.D-i  cxpOk.x  +  jkyy)dkx  dky 

-  M  .M 

(5o) 

mm  m 

~W  =  if f  k«  F*«  CXP^«  *  *  ikyy>  dk*  dky 

(5p) 

an; 

=-jk  sin#  T||  exp(-jkx  sinfl)  |  1  +  jk  z  cos*) 

(5q) 

an: 


dy 


KGri  r4'A'  +ikyy>dk»dky 

kyDX|  exp(jkxx+jkyy)dkxdky 


(5r) 


Substituting  equations  (5)  into  the  boundary  conditions  of  equations  (4)  and  the  two  di 
vergence  equations  and  evaluating  the  result  at  •*  =  s(x,y)  yields  the  following  six  alge¬ 
braic  equations  with  six  unknowns: 


(6a) 

=  <*, +  q*  ' 

(6b) 

\Kt  =  MVU*.. 

(6c) 

(6d) 
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Vkk'1  +  T|ky k'2Fxl  +  i}kfk'2Gyl  -  tjkaky Ft;  +  nk2k;G;,  +  Wkk'J  =  0  <6e) 

-tjk‘a kE Dm1  -  fjk’k, Ffl  -  t*2^,  +  tjkJkx Ft',  +  kk'2S0  =  0  <6f) 


where:  j(k2  -  k  * )  ( 1  -  R,.)  cos#  S(k  +  ksin0,  k  ) 

<>>= - 4^k - - 

V  =  -  jT||TjksinJff(ka  -  k'2 )  S(kx  +  ksinfl,  ky  )/2>rk'2 
W  =j(kx  +  k  sin0)  j?T||8infl  (k2- k’J)  S(kx  +  ksin$,  ky  )/2wk'2 


S0  -jkyt?T||  ginfl  (k2  -  k  2)  S  (kx  +  ksinfl,  ky  )/2irk2 


S(kx  +  ksinfl,  ky)  = 


g(x,y)  exp  [-j(kx  +  kgin0)x-jkyy]dy  dx 


f«j  ,j  , 

Use  hag  been  made  of  the  relation  (k  cos 0  -  R||  k  cos0)/k  =  Tj|k  cog 0  .  For  this  anal- 
ysis,  interest  is  only  in  the  fields  above  the  interface  so  that  solutions  need  only  be 
found  for  the  terms  Dxl ,  Gyl ,  and  FxJ  .  Using  equations  (6),  solutions  can  be  obtained 
for  these  three  terms: 


D*.  =  a 


Bo«„  - 


11  *22  "  a!2  a2l 


yi 


_  Ao  "m  -  Bo  a 

"  ai.  a22  -  a,2  a 


12^ 

21 


F_.  = 


■u 


k,G,+k.U.i 

IT - 


where 


\  =  -ukV^k.Q,  -tj^VQ,  -Wkik’ikk'2  -  Vk^kk*2 


B„  =  klk’lkk’2  S0  -tjk2  kxkykx  Q, 

a„  =  „l£k;k'2  +T7k;kk'2+^k2k2yki  +nk2k>x 
a,2  =  nkyk;k’2kx  +  i?  k2kyktkx 


a2i  =  nk’2kxk;ky  +t,k2kxkfky 

a„  =  nk'2^^  +  Tjk^k^k  +  r,k2kxk;2  +  „kV  kx 
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APPENDIX  B 


DERIVATION  OF  AVERAGES  USED  TO  CALCULATE 
MEAN  BACK  SCATTER  ED  POWER 

This  appendix  presents  the  derivation  of  the  necessary  averages  that  were  used  to  deter¬ 
mine  expressions  for  mean  powers.  The  two-dimensional  gaussian  distribution  of  two 
random  variables  z  and  z  with  zero  means  is 


C  is  the  autocorrelation  coefficient; 
a1  is  the  variance. 


z  -  2Czz'+  z 
2oJ(I  -C2) 


The  two-dimensional  characteristic  function  X  (v, ,  v2 )  of  the  distribution  p(z,z  )  is 
X  (v,  i  v2  )  =  ^ eJV|  1  *  p  (z,z)  dz  dz' 

-  oo  -oo 

X  (V, ,  v2 )  =  exp  O  (V;  +2C  v,  Vj  +  vj 
The  first  function  that  will  be  averaged  is 

ejki(,-,’) 

Where  k,  is  an  arbitrary  function  of  the  angle  of  incidence  0 

oo  oo 

<cik,(I  *^>  =  J* J"  e,k‘  Ie",k|  *  p(z,z’)dzdz 

- oo  -oo 

By  analogy  with  the  characteristic  function  x(v,  ,v2 ),  it  can  be  seen  that  v(  -  k,  and 
v2  -  k(  .  The  required  average  is  then 

<e,ki(*-0>  =  eXp|.  k|o2(l -C)| 

|l  ^ 

In  order  to  calculate  <Zxe  1  ’  '>,  it  is  necessary  to  rewrite  the  quantity  inside  the 

brackets  as 
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<Z,«*.«-C»=jjL  £  <«p{jk,(Z-Z')}  > 

=  {«pHs°’<i-c)]}|£ 

where  u  =  x  -  x’  and  v  =  y  -  y 

<Z>e*,<z-2  >>  =  .ik,o’  |£  expl-kjo’d  -C)|  . 

The  average  <Zye  *'  *  ’>  can  be  computed  in  exactly  the  same  manner  as  above, 
but  the  work  is  done  with  y  and  v  instead  of  x  and  u.  When  this  is  done,  the  result 
becomes 


<ZyeMt«(I'1)"- 


>-  ' iki°a  I7  exP  l-^0’  0  -C>1  • 


To  calculate  the  average  <ZX  Zx^k,<1'1  )>,  it  is  necessary  to  rewrite  the  quantity  in  the 
brackets  as 


jlij  3x  ]EJ  3x1 


<Z  Z'e,k'(**l  )> 


d  l,i  1  JC  ,  ,  j  7  . , 

du  ik,'  aiT  '•’■P  I- k.  “  <• 


-Ql} 

®*<l  -C)} 


The  rest  of  the  averages  which  depend  solely  on  the  large  surface  undulations  can  be 
computed  in  a  manner  similar  to  the  above  methods  and  the  results  will  be  stated  here 


<Z  Zejk'<Z'Z)>  =  <Z  Z  eik'(Z'Z)>  = 

x  y  >  y 


exp  |-k’oJ(l  -C)J 
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<z;«^z-rf>  =  <z<e*l<!1,>  =-jk,  o’|£exp|-k|c’<l-C)| 
=<^,<I-n>  -jk,»’  |£«p[.k1V<l-C)| 

<Z,Z;  *  k,’a'(|£)  J«xp[-k,’<-’(l-C)| 


The  averages  that  occur  with  the  slightly  rough  surface  s(x,y)  are  much  more  complicated, 
however,  if  we  could  compute  <SS#exp(jas  -  ja*s’)>  ,  all  the  other  averages  can  be  deter¬ 
mined  by  methods  used  for  the  large  undulations  alone.  The  Fourier  transform  of  s(x,y) 
is  S(kx ,  ky )  and  can  be  written  as 

OO  Oft 

S(kx,ky)  =  ^  J  J  s(x,y)e*ik**  e"V  dydx 

-  Oft  -  M 

<SS*exp(jas  -  ja*s’)>  =  <-^t  JJ ss'e",kxX  e^*  e  JkyV  e*** y  dydxdy'dx' 

•  exp(jas-  ja*s)> 

=  <^TT J  JJ j’  s.s^'e*’  e'ikyy  e^kyy  exp(jas3  -  ja*s4)dydxdy'dx> 
s,  =  s(x,y)  Sj  =  s(x',y)  &,  =  s(jf ,  y")  s4  =  s(x”',  y'") 

-  “jr  ^ J" ^ J"  e',k**  eik**  e>kyV  i,Jkyy  <s(  s2exp(ja83  -ja*s4)>  dydxdy’dx' 
The  multivariate  gaussian  characteristic  function  is  given  by  the  following  equation: 

M (jv, .  jv2 x  jvj ,  jv4 )  =JJJJ  eVi  **7*1  4  ^3*3  4  N*4  p(S(  ,sJf^,84)ds|  ds2  ds3ds4 

the  quantity  pfSpSj,  Sj,  s4)  is  the  multivariate  gaussian  distribution  function.  From  Daven¬ 
port  and  Root, ,J  the  multivariate  characteristic  function  can  also  be  written  as  the  follow¬ 
ing  equation  when  the  means  of  all  four  variables  are  zero: 


M(jv,,jv2,jv3,jv4)  =  exp 


(-  S  XmnVnVm) 

^  m  *  I  n  *  I  J 


\  -00 
mn  n  m  n 


’Davenport  and  Root,  An  Introduction  to  the  Theory  of  Random  Signals  and  Noite,  McGraw-Hill,  1958,  p.  153. 


85 


a^r  o 


=  <8,  \exp(ivaS  +jv4s4)> 

^,'(jV1’ jvj«  jv3,  jv4)  =  exp  {-  4*X.1V.2  +  V2V.  +X.3V3V.  +  X,4V,V4 

+  V.V2  +  Xa3V3+X»V4+V.V3+XMVa 

+  X3,V3  +  X34V3V4  +  X4,  V4  +  *41  V,  V,  +  X^V,  +  X„  vj  |J 

When  m  =  n,  =  1.  For  a  stationary  process  a*  =  o2  =  o*  =  •  Then  X,,  = 

\i  =  X3J  =  X44  1 

M(jvr  jv2,  jv3,  jv4)  =  exp^X,,  (vf  +  v,J  +  vj  +  v4V  2X,2  v,  v, 

+  2X13  v,  v3  +  2Xm  v,  v4  +  2XJ3  v2  Vj  +  2Xj4  v2  v4  +  2XJ4  v,  v4  |  J 

2 

a*  3v,  =  ‘  xn  exP  {"  4*1  xn  (v!  +  v!  +  VJ  +  \  ) +  2x„  v,  v2  + 

2X,3V,S  +  2x.4V.V4  +  2Xj j  v2  Vj  +  2Xj4  v,  v4  +  2XJ4  Vj  v4  j  | 

+  IV*X'<  +  XH  Vj  +  XI3  V3  +  X14  V4  I  (  Xn  Vj  +  Xu  V,  +  Xjj  Vj  +  XJ4  V4  | 

•  eXP  {  -tIX«  +  V3  +  s  +  V4)  +  2X,J  V,  Vj  +  2X,J  V,  Vj  + 

2Xm  v’i  v4  +  2Xj,Vj  V,  +  2Xj4VjV4  +  2XJ.VJV,  )  J 

dv  dv  {  +  XI2  ^Xl  3  V3  +  X14V4)<X3JV3  +  ^4  V4  )  1 

1  1  v|  =  v2  ”0  1 

'  cxp{‘ilx,.(S  +V4,)  +  2XJ4v,V4)| 


I»I  =V,=0 


Vj  =  a 


v4  =  -  a’ 


XIJ  °l  ^'IJ  ’ 

XI3  °l  ^IJ  ’ 

XI4  °l  ^14 

-  al  ’ 

^4  “  °l  ^4  * 

XM  ~  °l  ^34 
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?>■ 


It  can  easily  be  seen  that  the  term  (X„  ^  +  Xl4  v4  XX,,  v,  +  X,4  v4  )  is  of  the  order  of  cj* , 
while  the  term  X„  is  of  the  order  of  o,  .  The  variance  of  the  slightly  rough  surface 
(o,1)  will  always  be  smaller  than  one  for  the  wavelengths  to  be  considered  in  this  report. 
Therefore,  the  term  (X„v,  +  X|4v4)(X,,v,  +  X,4v4)  can  be  neglected  with  respect  to 
X„  .  The  above  equation  can  then  be  written  as 


avjdv, 


v  *  v  *  0 
1  v2  w 


_  =  °i  C12exp  -j-  [o*  (aJ  +  a*1)  -  2aa*C34ff*)|' 


<SS#  exp(jas-  ja#s)>  »  -—jr j  Jjf  e  ***  *****  ****  -  x’.y*/) 

•  dydxdy'dx'exp|-  -^o^fa1  +  a*1  -  288*0^11 
x-x'=u  y-y'=v  C,j  (u,  v)  =  C,  (u,v) 

<SS#exp(jas-ja*s)>  •  e*ik*ue’Jk>'v  e,*’(li'k»)  eiy(lV’k»)  oj* 

•  C(  (u^Jdudvdx  dy1  •  exp  |  1  a1  +  a*1- 288*0^  1 1 

since  the  surface  is  the  same  in  all  cases  with  respect  to  correlation  Cj,,  =  C,  (u,v) .  Inte¬ 
grating  the  above  equation  in  x’  and  y\  the  Dirac  delta  functions  appear. 

<SS*exp  (jas  -  ja*s)> 

m  m 

6(k|(  -kx)S  (k^  -  ky)a, *  J  J  C,  (u,v)  e'ik*u  e'**"  dudv 

•  exp  | ^ o,5  [aJ  +  a*J  -  2aa*Cj| 

The  double  integral  in  the  last  equation  represents  the  roughness  spectrum  W  (kx  ,ky): 

^(kx,ky)=  ^  J  J  C,  (u,v)c',k*u  e'JV  dudv 

-M 

The  final  average  is  then 

<SS#exp(jas  -  ja#s)>  *  2s6  (k^  -  kx )  6  (k^  -  ky )  o*  W  (kx  ,ky ) 

•  exp  y  o’  |aJ  +  a#J-  2aa*r.(  |  j 
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I 


When  C,  (u,v)  is  gaussian  and  of  the  form 
C,  (u,v)  =  e-("1+^'a 


then  W(kx,ky)i« 


W(k* ’M ‘  T«p{-  fk*  +  ky  JP/4J 

klS  derjV<,ti0n8  «  **  analy8i8  8eCti°"  °f  the  rePort,  k  must  be  replaced  by  k  + 

ksmd  and  is  actually  approximated  with  ky +k  sin  0  .  P  Y  * 

,mtst  <,-ss*'xp('*-  ^  -  ■» *•»«■*-* 
<sxSS*cxp(jas-ja*8'>>  * ^  <SS*exp(jas  -  ja*s')> 

*  ]T  ■^J’<SS*exp(ja8  -  ja*8'>>  -|^- 

<sx  SS*exp  (jas  -  ja*s')>  *  -  ±  {  2*6  -  ky)  6  (k  -  kx  ) 

•  W(kx,ky)oIiexp|-  i  o/fa2  +a#1  -2aa*C,  J  J 


.•iS 


ja#  °*  *T  2ff6  <k'r  k* > 5  (k’y  -  ky )  w  (kx ,  ky ) 


■  CXP  f  y  (a2  +  a*1  -  2aa*c,  J  j 

Using  this  same  technique,  the  following  averages  can  be  computed  easily: 


<sySS*exp(jas-ja*s')>  *  -j2*aViS. 


y  w°'Jd  ) fi  (ky  - ky) 

•  w  (kx  <  ky )  exp  i-02  fa2  +  a*2-  2aa*C,  ]} 
<8*^S*  exp  (jas  -  ja*s')>  *  ~  j2rra  o*  ^6  (k  -  k, ) Mk  -  ky ) 

'  W(k„,ky)exp|-  I0|2  |a2  +  a*2  -  2aa*C,  )J 
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4  3C. 

<8;SS*exp(jas-jaV)>  *  -j2irao1  fi  (k;  -  k,)  6  (k'y  -  ky ) 


•  w(kx .  ky)  e*P  \  o,1  [■*  +  «#I  -  2aa*C, )  J. 
<sx  8^  SS#exp  (jas  -  ja#«)>  »  -  2no\  S  (k'^  -  k% )  6  (k^  -  ky  ) 

■  **.■>.>  {0  •  s?f 

•  exP  { -  °!  (*5  +  a*J  -  2aa#C,  1 1 

<8y8'ySS*exp(jas-ja*8)>  **  -  2aoJ  6  (ky  -  ks)  6  (k'y  -  ky  )  W  (kK  ,ky ) 

■  +  “■  “•  (t^)} exp  {  '  T  |a’ +a*’ 

-  2aa*C(  1 1 

<sxsy  SS*  exp  (jas  -  ja*s)>  =  <sy  SS*exp(jas  -  ja*s)>  -  2jto* 

■  »<k;  -  k.) «  (k,  -k,)  W  (k,.  k,)|~  +  (^) 

•  exp  -j  o*  (a*  +  a*’  -  2aa*C,  J  | 

All  the  required  averages  needed  to  derive  average  power  have  now  been  determined. 


APPENDIX  C 


EVALUATION  OF  INTEGRALS 


This  appendix  shows  the  derivation  of  some  of  the  integrals  needed  in  the  analysis 
section  of  the  report.  The  two  basic  parameters  which  are  used  in  the  needed  integrals 
are  a  and  r.  There  is  only  one  integral  in  r  that  needs  to  be  specified  and  this  can  be 
obtained  from  Gradshteyn  and  Ryzhik. 14 


One  of  the  integrals  in  a  which  must  be  calculated  is: 


/in 

cos  (or)  e*^,in  ( 


Let  u  =  j0  sin  a 

du  =  j(3  cos  a  da 


•j0tina  ,  /  -u  du  1  /  -u  .  n 

'  ,la  J 


cos (a) e 


Using  the  above  result  will  permit  the  evaluation  of  the  integral: 


/in 

-J0iin  ( 

sin  a  e 


a  -  -je  +jcosa 


sin  a  i 


-10 fin  a 


da  =  -j  J"  c'V^da  +  j  J 


-jflilna  , 

cos  a  e  da 


r” 

J  sin  a  e‘^*ina  da  --j2rj  ( 


1.  S.  (irifWiteyn,  I.  M.  Kyzhtk,  Tablet  of  Integrals  Series  and  Products ,  Academic  Preaa,  Inc.,  1965. 


vo 


J_,  (-0  =  -J,(-«  “J.W 


The  J  functions  are  Bessel  functions  of  the  first  kind: 

/in 

sin  a  e*^,ln1do  =  -j2irj1(0) 
o 

The  next  integral  which  must  be  evaluated  is 

/in 

i  -jfliln  a  , 
cos  a  e  r  da  . 

o 

This  integral  can  be  worked  by  using  the  method  of  integration  by  parts: 


u  -  cos  a 


du  =  -  sin  a  da 


dv  =  cos  a  e'J^,lnada 
'  in 


v  =  }c 


-j0«in  a 


/s  -j0»ina  .  _  icosa  -J0«ina  2n  ,  \  f 

cos  ae'K  da  =  -  e  +  J  I 

®  oo 

m  r  ** 

j“  cos,a  c  ^,ina  da  =  0  +  -jj- J 

0  o 

/Jff 

I  -J0iina  .  Jj  (0) 

cos  a  r  da  =  - ! - 

a 


in 


sin  a  e 


-J0>in  a 


-  j/3im  a  , 
sin  a  p  da 


Another  integral  which  needs  to  be  evaluated: 


/in 

sin’ 


10  tin  a  , 

ae  da 


/air  *>in  fin 

■  1  -J0«in  a  J  /  -)0un  Q  ,  /  3  -)0 

sin  a  e  da  =  I  e  da  -  f  cos  a  e 


-jpiin  a 


da 


=  2*  IJ„W  •  ^ 


=  2.  [M*>.J,W] 


da 
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The  last  integral  which  needs  to  be  evaluated  is: 


cos  a  sin  a  e 


f,<Utaada 


This  integral  can  be  solved  by  using  the  method  of  integration  by  parts: 


u  =  sin  a 


j  _  -J0sln  a  i 

dv  =  cos  at  da 


du  =  cos  a  da 


v  =  j  « 

V  Je 


sin  a  cos  a  e 


/lit 

- 1/3  tin  a  , 

cos  a  e  K  da 


-J0iln  a  ■  n 
sin  a  cos  a  e  da  =  0 
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